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Chapter  1 


RADIAL  PRESHEATH  MODEL  FOR 
FLOWING  PLASMAS 

1.1  INTRODUCTION 

The  transverse  interaction  of  a  collisionless  plasma  with  two  dielectric  walls  confining  it  was 
treated  originally  by  Tonks  and  Langmuir[54],  who  solved  the  fundamental  issues.  The  zero- 
current  condition  at  the  walls  means  that  fluxes  of  ions  and  electrons  to  the  walls  must  be 
equal,  which  implies  a  potential  profile  decreasing  monotonically  from  around  the  center  of  the 
channel  towards  the  two  walls.  Tonks  and  Langmuir  showed  that,  in  the  zero  Debye-length 
limit,  the  plasma  structure  consists  of  two  thin  non-neutral  sheaths  tied  to  the  walls  and  a 
quasineutral  region  (the  presheath)  occupying  the  bulk  of  the  channel.  Each  type  of  region  is 
solved  separately  in  their  own  distinguished  scale.  The  condition  for  a  correct  presheath/sheath 
transition  is  unique  and  consists  in  the  plasma  flow  being  sonic  there  (the  Bohm  condition [20]). 
Apart  from  the  three  basic  wall  geometries,  Tonks  and  Langmuir  (interested  in  thin  cylindrical 
probes  immersed  in  cylindrical  channels)  discussed  the  case  of  an  annular  geometry,  pointing 
out  the  lack  of  symmetry  between  the  profiles  of  the  inner  and  outer  presheaths  (these  ones 
separated  by  the  point  of  maximum  electric  potential). 

A  central  contribution  of  Tonks  and  Langmuir  is  the  plasma  balance  equation ,  which  estab¬ 
lishes  the  plasma  production  rate  required  to  sustain  a  steady-state  discharge.  This  production 
rate  is  proportional  to  the  channel  width,  as  it  comes  out  from  equating  the  volumetric  produc¬ 
tion  of  plasma  to  the  losses  at  the  walls.  In  Tonks-Langmuir  model,  where  plasma  production 
is  dominated  by  ionization,  the  plasma  balance  equation  fixes  the  ionization  rate  and,  as  a 
consequence,  the  plasma  temperature;  indeed,  a  steady-state  discharge  is  not  possible  when 
the  plasma  production  rate  is  larger  than  the  maximum  ionization  rate  of  the  gas.  Tonks  and 
Langmuir  considered  different  radial  profiles  for  the  ionization  rate,  but  it  was  shown  later  that 
the  mean  value  of  the  plasma  production  rate  is  independent  of  its  radial  profile  [33,  23]. 

Bissell  et  al.  [16]  (reviewing  several  other  papers)  and  Scheuer  and  Emmert[50]  propose 
macroscopic  formulations  of  the  Tonks-Langmuir  problem  (for  planar  geometry  only),  with 
different  closure  hypotheses  for  the  ion  energy  equation.  They  conclude  that  macroscopic 
models  can  approximate  well  the  kinetic  results,  mainly  in  relation  to  particle  and  energy 
fluxes  towards  the  walls. 

As  the  Tonks-Langmuir  model,  these  fluid  models  include  ionization  as  the  only  source  term 
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in  the  plasma  equations  and,  therefore,  are  restricted  to  no-flowing  plasmas  basically.  When  the 
plasma  flows  along  the  confining  walls,  the  problem  becomes  two-dimensional(2D),  with  mutual 
interaction  between  the  longitudinal  and  transverse  (i.e.  axial  and  radial)  dynamics.  However, 
at  each  axial  cross-section  (of  a  long  channel)  we  may  consider  the  one-dimensional(lD)  problem 
that  describes  the  radial  dynamics  of  the  plasma.  The  requires  to  take  into  account  the  influence 
of  the  axial  flow  on  the  radial  motion.  In  the  macroscopic  formulation  this  can  be  modeled 
with  adequate  source  terms  in  the  radial  equations. 

The  goal  of  this  paper  is  to  generalize  the  presheath/sheath  model  of  Tonks-Langmuir 
to  plasmas  flowing  axially,  with  emphasis  on  the  coupled  influence  of  ionization  and  axial 
flow  on  the  radial  solution.  In  a  preliminary  study  (with  planar  geometry  and  cold  ions)  we 
showed  already  that  losses  to  the  wall  are  reduced  by  an  effective  radial  friction  due  to  axial 
variations  of  the  radial  conditions[9].  Here  we  develop  the  general  model  for  warm  ions  and 
annular  geometries.  Thus,  planar,  cylindrical,  and  annular  channels  are  covered;  the  differences 
between  the  inner  and  outer  presheaths  in  the  latter  case  are  an  extra  aspect  of  interest.  No 
secondary  emission  at  the  walls  will  be  considered  here,  so  that  the  sheath  problem  reduces  to 
a  well-known  solution  and  the  attention  of  the  paper  is  focused  on  the  presheath  model. 

An  example  of  application  of  the  present  work  is  the  plasma  discharge  in  an  annular  Hall 
thruster.  Since  an  accurate  two-dimensional  analysis  of  this  problem  is  unavailable  yet,  partial 
analyses  of  the  radial  and  axial  problems  are  of  great  interest.  In  particular,  the  knowledge  of  the 
radial  plasma  structure  would  allow  one  to  estimate  the  losses  at  the  lateral  walls  at  each  cross- 
section,  which  could  then  be  included  as  sink  terms  in  ID  axial  models  of  the  discharge[26,  11]. 
The  2D  Hall  thruster  discharge  has  been  chosen  here  as  the  basis  to  formulate  the  ID  radial 
model  for  a  flowing  plasma.  However,  we  believe  this  formulation  to  be  general  enough  to  apply 
to  other  flowing  plasmas  interacting  with  two  walls.  In  order  to  emphasize  this  last  point  and 
to  free  the  radial  model  from  the  Hall  thruster  case,  the  derivation  of  the  ID  radial  model  from 
the  2D  axisymmetric  model  of  the  discharge  has  been  confined  to  Appendix  A. 

In  a  recent  paper  on  a  2D  model  of  the  Hall  thruster,  Keidar  et  al. use  a  presheath/sheath 
model  for  the  lateral  walls[37].  They  claim  that  for  a  smooth  presheath/sheath  solution  the 
radial  ion  velocity  at  the  sheath  transition  must  be  smaller  than  the  (sonic)  Bohm  velocity,  in 
contradiction  to  the  well  established  sheath  theory  between  a  wall  and  a  plasma[54,  48].  The 
main  arguments  supporting  the  sonic  transition  are  recalled  in  Appendix  B  here.  Comments 
on  other  aspects  of  Keidar’s  paper,  like  their  results  on  the  interaction  between  the  axial  and 
radial  responses  are  out  of  the  scope  of  this  paper. 

The  rest  of  the  paper  is  organized  as  follows.  In  Sec.  II  we  formulate  the  radial  model,  the 
boundary  and  compatibility  conditions,  and  the  integration  procedure.  In  Sec.  Ill  we  discuss 
solutions  for  different  cases  and  geometries.  Conclusions  are  presented  in  Sec.  IV. 

1.2  MODEL  FORMULATION 

1.2.1  The  quasineutral  presheath 

Figure  1.1  sketches  the  situation  on  an  axial  cross-section  of  an  annular  chamber.  The  plasma 
is  supposed  to  flow  axially  between  two  long  dielectric  walls  placed,  in  cylindrical  coordinates 
(x,r,6),  at  r  —  rw1  and  r  —  rw  —  T'w  +  h,  with  h,  the  channel  width.  We  consider  the  zero 
Debye-length  limit,  Xd/h  — »  0,  allowing  a  two-scale,  presheath/sheath  analysis,  where  h  and  A d 
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are  the  distinguished  scales  of  the  presheath  and  sheath,  respectively.  From  Appendix  A,  the 
pertinent  ID  fluid  equations  in  the  presheath  are 


1  d 

7Z~'TTliVri  — 

r  or 

(i.i) 

dvri  _  d  d(p 

miTliVri  —  [Tlilij  CTli 

or  or  or 

(1,2) 

\  c)  3  T'  n  0  1 

(r  TniVri)  +  (rvri)  =  m mvI  +  vtn{Th 

r  or  2  r  or  2 

(1,3) 

ne  —  const  x  exp (e</>/Te), 

(1,4) 

ne  =  Hi. 

(1,5) 

Here:  I'j,  uw,  vr.  vt  are  frequencies  for  ionization,  ion  net  production,  effective  radial  friction, 
and  axial  energy  contribution,  respectively,  defined  in  Eq.  (A. 7).  We  assume  these  frequencies 
to  be  known  functions  of  r  in  any  axial  position.  Moreover,  we  restrict  the  discussion  here  to 
Vj{r)  —  const,  (j  —  i,w,r,t ),  because  there  are  no  fundamental  novelties  when  treating  more 
general  profiles  (see  Ref.  [23]  and  comments  in  Sec  III. A). 

Operating  with  the  above  equations,  the  radial  problem  in  the  presheath  reduces  to  integrate 
the  two  coupled  equations 


(Te  +  | Ti  -  miV^  =  (vr  +  j)miV ft 

-  (Te  +  ^T)  —  +  vw(T  +  1^)  +  (^  -  \vw)Th 
dr  3  3  3 

3  dT  1  2  rp  ( dvri  vri  \  .  5  . 

2Vrilh  =  Ui2miVri  ~  Tl \lfr  +  V )  +  VwTl  +  ^  "  2Vw’Th 

plus  one  differential  equation  for  the  electric  potential, 

d  e<p  Id 

h.Trb  =vw- 
or  I e  r  or 

which  corresponds  to  the  so-called  plasma  equation  in  kinetic  models[54]. 
Equation  (1.6)  shows  that  radial  derivatives  become  singular  at 


Vri  —  Cs  (Ti ) 


lTe  +  lTt 


3  * 


rrii 


(1.6) 

(1.7) 


(1.8) 


(1.9) 


with  cs  the  local  sound  velocity  of  the  plasma.  In  particular,  the  electric  field,  d(j)/dr ,  becomes 
infinite  there,  indicating  the  failure  of  the  quasineutral  approximation  and  the  transition  from 
the  large  h-scale  to  the  thin  A^-scale  of  the  sheath.  Therefore, 


VriQ  —  CsQt  V riQ '  —  (‘sQ1  (1.10) 

are  the  Bohm  conditions  defining  the  transition  points  (Q  and  Q'  Fig.  1.1)  to  the  outer  and 
inner  non-neutral  sheaths.  In  Appendix  B  we  show  that  these  transition  conditions  are  the 
only  valid  ones. 
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The  electric  potential  must  present  a  maximum  at  some  intermediate  position  of  the  channel 
(point  M  in  Fig.  1.1),  d<f>/dr\M  =  0.  A  local  expansion  of  Eqs.  (1.6)-(1.8)  around  point  M  yields 


vri(r)  =  vw (r  -  rM)  +  0[(r  -  rM )2], 


Ti(r )  =  Tm  +  0[(r  -  rM)2], 


(1.11) 


and  the  compatibility  condition 

vt  =  5vw/2.  (1.12) 

As  Eq.  (A. 9)  in  Appendix  A  confirms,  this  condition  corresponds  to  the  correct  balance  between 
the  radial  and  axial  energy  fluxes  for  ions. 

Condition  (1.12)  makes  identically  zero  the  last  term  of  each  of  the  equations  (1.6)  and 
(1.7).  Then,  for  the  zero-ionization  case  —  0),  Eqs.  (1.1)  and  (1.3)  yield  the  isentropic  law 


T 

1 1 


—  const, 


(1.13) 


which  substitutes  to  Eq.  (1.7). 

Convenient  dimensionless  variables  to  integrate  Eqs.  (1.6)-(1.8)  are 


r 


-  T-  =  —  6=  6 —  —  ^ 

h*  i  Te  *  Te 


Vri  — 


VTe/mi 


V3  = 


Vj  ,  . 

—  U  =  hw,r,t), 
Vo 


(1.14) 


with 

Vo  =  h~X\jTe/rrii  (1.15) 

the  characteristic  frequency  for  radial  processes.  The  convenient  integration  procedure  departs 
from  point  M  towards  points  Q  and  Q',  using  the  expansions  of  Eq.  (1.11)  as  initial  conditions. 
The  two  independent  integrations  end  where  Bohm  conditions  (1.10)  are  satisfied,  and  deter¬ 
mine  fq  and  tq>  in  terms  of  five  input  parameters:  T)m  ,  rM,  vr.  and  i>w,  while  Dt  satisfies 
Eq.  (1.12). 

The  sheaths  thicknesses  are  of  the  order  of  the  Debye  length.  Thus,  in  the  zero  Debye-length 
limit,  the  sheaths  are  seen  as  two  surface  discontinuities  in  the  presheath  scale: 

Tw-rq,  rw>  -  rQ> .  (1.16) 

Since  the  channel  width  is  known,  the  presheath  solution  must  satisfy  the  condition 


?q  ~  fQ>  =  1. 

This  implies  an  extra  relation  among  the  above  five  input  parameters.  After  exchanging  rM  by 
the  annular  ratio 

/3  —  r  w'/rwi 

as  known  parameter,  the  extra  relation  leads  to  a  relationship 

Vw  vw  (yi ,  vr ,  Tjjvf,  /3) .  (1.17) 
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This  is  the  generalized  form  of  the  Tonks-Langmuir  plasma  balance  equation  for  flowing  plasmas 
and  annular  geometries;  the  limits  ft  — >  1  and  fi  =  0  correspond  to  the  planar  and  cylindrical 
cases,  respectively. 

The  comparison  of  the  present  model,  Eqs.  (1.1)-(1.5),  with  those  of  Refs.  [16]  and  [50],  shows 
that  these  ones  cover  only  the  cases  i)  /3  —  1  (i.e. ,  a  planar  geometry)  and  ii)  vr  —  ^  —  vw  (i.e., 
a  no  flowing  plasma).  Therefore,  the  four  free  parameters  of  Eq.  (1.17)  are  reduced  in  their 
case  to  just  one:  Also,  these  no-flowing  models  use  z^|Tn  in  Eq.  (1.3)  instead  of  z^T*, 

which  agrees  with  Eq.  (A. 7)  for  the  no-flowing  case.  This  difference  in  the  heating  term  has 
interesting  consequences  on  the  balance/compatibility  condition  for  the  internal  energy.  In  our 
model,  the  term  z^T*  leads  to  the  compatibility  condition  (1.12),  which  states  that  the  radial 
heating  frequency,  ist,  depends  on  the  radial  production/loss  frequency,  uw.  whereas  the  value 
of  TiM  is  free  (which  agrees  with  the  idea  of  depending  partially  on  the  axial  discharge 
conditions).  On  the  contrary,  in  the  no-flowing  models  of  Refs.  [16]  and  [50],  which  do  not 
include  the  heating  frequency  vt  (that  is,  they  neglect  any  axial  contribution  to  ion  heating), 
the  balance/compatibility  condition  fixes  the  value  of  TiM  (as  proportional  to  Tn).  [Differences 
in  the  definitions  of  the  ion  temperature  in  the  different  models  have  no  relevant  consequences 
and  have  then  been  ignored  in  the  preceding  discussion.] 

1.2.2  The  sheaths 

The  role  of  the  two  radial  sheaths,  W'Q  and  QW  in  Fig.  1.1,  is  to  provide  the  potential  drop 
needed  to  maintain  the  zero-current  condition  at  the  wall.  For  the  zero  Debye-length  limit,  the 
sheaths  are  collisionless  boundary  layers  perpendicular  to  the  walls.  The  2D  Poisson  equation, 
Eq.  (A. 11),  reduces  to  the  ID  radial  form,  Eq.  (B.l),  and  plasma  equations  (1 . 1 )- (1 .5)  become 
algebraic,  conservation  relations  for  the  fluxes  of  different  plasma  magnitudes.  In  particular 
and  taking  the  case  of  the  outer  sheath,  the  ion  density  flux  to  the  outer  wall,  griw  =  (^Rn)w, 
satisfies 

9riW  —  9riQ  —  ^iQ^sQi  (1.18) 

with  riiQ  =  riiM  exp(— e<f>QM/Te)  and  4>qm  =  4>m  ~  4>q  the  potential  drop  at  the  presheath. 

For  no  secondary  emission  at  the  walls  and  a  semi-Maxwellian  distribution  of  the  electrons, 
the  zero  current  condition  follows 


(  efiwQ  \  [ 

9riQ  =  neQ  exP  ( - 


27 rme  ’ 

with  (fwQ  —  4>q  —  the  potential  drop  at  the  sheath.  Solving  for  4>wq ,  one  has 


(1.19) 


e<t>WQ  _  1  j]n 


771; 


TP 


27T7TZ, 


-Mi  +  lff)]. 


(1.20) 


For  all  practical  values  of  Tjg/Te,  the  last  term  on  the  right  can  be  disregarded,  so  that  the 
sheath  potential  drop  is  practically  independent  of  the  presheath  solution. 

The  above  expressions  relate  the  main  sheath  parameters.  The  Poisson  equation  (B.l)  needs 
to  be  solved  only  to  determine  the  spatial  structure  of  the  sheath  and  to  find,  in  a  way  equivalent 
to  that  used  in  Appendix  B,  that  the  (outer)  sheath  solution  exists  only  for  vriQ  >  csq. 
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1.3  PLANAR  AND  ANNULAR  SOLUTIONS 


1.3.1  Planar  geometry:  cold  ions 

We  consider  first  the  simple  case 

P  — >  1)  TiM  =  0,  Vi  — >  0.  (1-21) 

The  first  condition  implies  that  the  channel  is  quasiplanar  (satisfying  1/r  <C  d/dr),  so  that 
the  presheath  is  composed  by  two  symmetric  regions,  with  =  (tq  +  tqi)/ 2  the  channel 
mid-point,  and  we  just  need  to  solve  the  region  MQ.  The  second  and  third  conditions  yield 
Tj(r)  =  const  =  0  as  the  solution  of  Eq.  (1.7). 

For  the  basic  ’no-friction’  case,  vr  —  0,  the  velocity  and  potential  profiles  [Eqs.  (1.6)  and 
(1.8)]  satisfy 

vri  -  Vri/3  =  {f  -  rM)vWi  <f>  =  -v*i/ 2.  (1.22) 

Setting  vriQ  —  1  at  ?q  —  ?m  +  1/2  [Eq.  (1.10)]  one  has 

vw  —  4/3,  (1.23) 


which  is  the  plasma  balance  condition  for  this  basic  case.  It  states  that  the  total  production 
frequency  is  fixed  by  the  channel  width  and  the  plasma  temperature.  Notice  that  the  plasma 
balance  condition  is  the  consequence  of  two  facts:  (a)  wall  electrical  properties  imply  the 
presence  of  a  thin  non-neutral  sheath;  and  (b)  the  transition  to  this  sheath  requires  the  ion 
flow  to  be  sonic. 

[Were  not  the  frequency  profile  vw(r)  constant,  (r— f^)^  in  Eq.  (1.22)  should  be  substituted 
by  JT  vw{f')df'  ,  and  the  balance  condition  would  yield  f~wf  vw{f')df'  —  4/3.  Since  vri  must 
be  zero  at  point  M  and  equal  to  ±cs  at  the  walls,  the  shape  of  vw(r)  modifies  only  the  local 
gradients  of  vri(r ):  when  vw  tends  to  concentrate  around  point  M,  vri  tends  to  present  an 
intermediate  inflexion  point[23].  In  any  case,  since  the  shape  of  vw(r)  depends  on  both  the 
axial  flux  and  the  ionization  rate,  very  steep  shapes  are  unlike.] 

For  the  positive  ’friction’  case  vr  >  0,  Eqs.  (1.6)  and  (1.8)  yield 


1  +  ar 


arctan  {vriy/a/) 


Vri  cxrvw(r  T m)  i 


1  +  cnr 
2  ar 


ln(l  +  OtrVri), 


(1.24) 

(1.25) 


with  ar  —  vr/vw  a  convenient  parameter [9].  Setting  vriq  —  1  at  f  =  rM  +  1/2,  one  obtains  the 
plasma  balance  condition  relating  now  vw  to  vr.  Figures  1.2(a)-1.2(c)  (case  T)m  =  0  and  Vi  —  0) 
show  the  evolution,  with  the  friction  frequency,  of  the  production  frequency,  the  potential  drop 
at  the  presheath,  and  the  plasma  current  to  the  wall.  For  vr  — >  0,  one  recovers  vw  =4/3  and 
4>qm  =  1/2;  for  vr  1,  one  has  vw  —  ti2 /vr  and  4>qm  —  In  vT.  These  results  indicate  that  larger 
ion  frictions  (i)  mean  larger  potential  and  density  drops  across  the  channel  in  order  to  bring 
the  ion  flow  to  the  sonic  condition,  vriQ  —  1,  at  the  sheath  tied  to  the  wall,  and  (ii)  yield  lower 
production  frequencies.  Therefore,  the  axial  motion  turns  out  to  be  an  efficient  mechanism  to 
inhibit  the  plasma  losses  to  the  wall. 
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However,  the  ion  axial  motion  can  lead  to  negative  values  of  Dr  too,  which  increase  the 
plasma  current  to  the  wall.  For  vr  <  0,  vri(f,  vr)  and  vw{vr)  satisfy 


Vri  T 


l+ar  1  —  nrj|o;r.|1/2 

2\ar\1!2  l+u„| a,]1/2 


ar\uw(r  ~  fM), 


(1.26) 


and  Eq.  (1.10).  These  relations  show  that  dvw/d(— vr)  >  0,  but  solutions  with  negative  friction 
exist  only  up  to  vr  —  —2  (yielding  vw  =  2).  Were  —  vr  >  2,  the  plasma  velocity  would  not  reach 
the  Bohm  velocity,  Eq.  (1.10),  and  a  steady  radial  solution  would  not  form. 


1.3.2  Planar  geometry:  general  solution 

To  solve  numerically  Eqs.  (1 .6)- (1 .8)  for  >  0,  T \m  >  0,  and  0  =  1,  we  use  the  spatial 
variable  £  =  vw{f  —  rM)  and  the  ratios  oq  =  Vj / vw  ( j  —  r,i).  The  integration  from  points 
M  to  Q  yields  the  relation  vw  =  2  £q(ar,  ck,,  T)m),  where  from  the  plasma  balance  condition 
vw  —  vw(yr.K%M)  is  obtained. 

The  comparison  of  the  solid  and  dash-dot  lines  in  Figs.  1.2(a)-1.2(d)  illustrates  the  effect 
of  warm  ions  on  zero-ionization  situations.  First,  since  the  plasma  sound  speed  increases  with 
TiM,  one  has  that  vw  and  wall  losses  increase  too.  Second,  for  —  0,  T,t (r)  decreases  towards 
the  walls  following  the  isentropic  law  (1.13)  and  yielding 


TiQ/TiM  =  exp(— 30qM/2)  <  1. 


Figure  1.3  and  the  dashed  line  of  Figs.  1.2(a)-1.2(d)  illustrate  the  influence  of  i \  for  different 
values  of  Vi/vr.  Although  ^  contributes  to  vr  [Eq.  (A. 7)],  for  the  present  discussion  it  is 
convenient  to  take  them  as  independent  parameters.  Then,  Figs.  1.2(a)  and  1.2(b)  indicate 
that  vT.  and  not  i z,,  determine  the  net  production  frequency.  The  main  effect  of  is  to  increase 
the  dispersion  of  ion  velocities,  that  is  the  ion  temperature.  Figure  1.4  shows  the  competition 
between  the  cooling  effect  of  ion  rarefaction,  Eq.  (1.13),  the  ion  heating  due  to  (and  increasing 
with)  ionization. 

The  no  axial- flow  case:  vw  —  vr  —  £j,  treated  by  Tonks  and  Langmuir  and  Refs.  [16]  and 
[50],  corresponds  to  a  single  solution  point  for  each  T \m]  these  are  the  asterisks  marked  in 
Fig.  1.3.  As  Tonks  and  Langmuir  explained  (for  TiM  —  0),  the  plasma  balance  condition  in  the 
no  axial- flow  case  determines  totally  the  ionization  frequency,  i and,  therefore,  the  electron 
temperature.  In  cases  where  the  axial  flow  is  relevant  (like  the  plasma  discharge  in  a  Hall 
thruster),  Te  is  determined  by  the  axial  energy  balance  mainly[10],  and  the  difference  between 
V{(Te)  and  the  total  production  frequency,  vw(Te,h)  is  provided  by  the  axial  variation  of  the 
ion  axial  flux  [Eq.  (A. 8)].  Finally,  notice  that  in  the  Tonks-Langmuir  problem,  a  steady-state 
plasma  discharge  requires  vw  to  be  less  than  the  absolute  maximum  of  ^(Te),  whereas  in  a 
plasma  flowing  axially,  the  steady-state  solution  relies  more  on  axial  conditions. 

1.3.3  Annular  geometry 

In  this  general  case  the  convenient  integration  variable  of  Eqs.  (1.6)- (1.8)  is  rj  —  fvw.  The 
separate  integrations  of  the  presheath  equations  from  point  M  towards  points  Q  and  Q',  which 
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includes  tjm  as  an  input  parameter,  yields  t]qi  and  t)q  as  functions  of  7]m,  cxr,  cq,  and  T)m-  Then, 
using  the  two  equalities 

Vw  =  Vq-  Vq>  =  Vq/?q,  (1.27) 

the  plasma  balance  condition  in  the  form  of  Eq.  (1.17)  is  obtained. 

Figures  1.5(a)-1.5(c)  depict  spatial  profiles  of  different  magnitudes  for  zero-ionization  cases 
(z Vi  —  0)  and  different  annular  ratios.  Point  M  tends  to  drift  towards  the  inner  wall  as  ft 
decreases,  yielding  two  asymmetric  presheaths,  with  steeper  gradients  in  the  inner  one,  as 
Tonks  and  Langmuir  deduced.  Plasma  gradients  tend  to  diverge  as  ft  — >  0,  but  the  solution 
cannot  be  continued  down  to  ft  =  0  (i.e.,  a  cylindrical  chamber),  since  then  point  M  coincides 
with  the  channel  axis  and  the  inner  presheath  disappears. 

For  ^  —  0,  the  ion  temperature  follows  the  isentropic  law  (1.13).  Consequences  of  this  and 
the  Bohm  conditions  (1.10),  are  vriQ  =  vriQi  and  riiQ  =  n^,  independently  of  ft.  Therefore, 

in  spite  of  the  presheath  asymmetry,  one  has  equal  wall  currents  and  potential  drops  in  both 

presheaths, 

Vi  —  0  :  griw  =  griWi  4>Q'm  —  4>qm- 

Figures  1.6(a)-1.6(c)  show,  for  different  ft,  the  evolution  with  v,t  of  the  main  parameters  of 

the  inner  and  outer  presheaths.  In  Fig.  1.6(a)  we  see  that  vw  is  affected  weakly  by  the  annular 
ratio,  ft;  notice  that,  similar  to  Fig.  1.3,  the  no  axial-motion  case  corresponds  to  just  one  point: 
vw  —  Vi,  for  each  /3-line.  Figures  1.6(b)  and  1.6(c)  show  that  for  ^  >  0  and  ft  <  1,  wall  currents 
and  potential  drops  are  different  in  the  two  presheaths.  Differences  increase  as  ft  decreases; 
departures  from  the  planar  solution  become  relevant  for  ft  <  0.7,  roughly.  Taking  as  a  reference 
the  planar  case,  potential  drops  tend  to  be  smaller  in  the  inner  presheath  and  larger  in  the  outer 
one  [Fig.  1.6(b)]  whereas  plasma  currents  [Fig.  1.6(c)]  follow  the  opposite  behavior.  The  sheath 
solution  (1.20)  indicated  that  presheath  asymmetries  of  an  annular  chamber  do  not  produce 
further  asymmetries  in  the  inner  and  outer  sheaths.  In  conclusion,  in  an  annular  chamber  with 
non-zero  ionization  the  particle  and  energy  fluxes  are  larger  at  the  inner  wall  than  at  the  outer 
one. 

Ion  temperatures  at  points  Q  and  Q'  for  the  cases  of  Fig.  1.6  are  about  0.4  —  0.5Te  and  their 
dependence  on  vr  and  ft  follow  the  trends  of  the  respective  potential  drops.  For  v,t  —  0,  vr  >  0, 
and  the  rest  of  conditions  the  same  as  in  Figs.  1.6,  results  are  very  similar  to  those  of  Figs.  1.6, 
except  for  the  ion  temperatures,  which  are  smaller  due  to  the  isentropic  behavior,  Eq.  (1.13). 


1.4  SUMMARY 

A  model  of  the  radial  structure  of  a  plasma  flowing  along  two  dielectric  walls  must  take  into 
account  the  contributions  of  the  axial  flow  into  the  radial  equations.  In  a  macroscopic  formu¬ 
lation,  this  can  be  done  by  including  three  terms,  representing  a  particle  source,  an  effective 
friction,  and  an  internal  energy  source  in  the  equations  of  continuity,  momentum,  and  energy, 
respectively.  The  production/loss  of  ions  may  come  from  either  ionization  or  the  axial  flow; 
the  effective  radial  friction  comes  from  both  ionization  and  the  axial  variation  of  the  radial 
velocity;  and  the  internal  energy  of  the  axial  flow  contributes  to  the  energy  balance.  Compared 
to  the  classical  model  of  Tonks-Langmuir  (and  its  macroscopic  versions)  the  present  model  has 
three  new  parameters  (or  degrees  of  freedom):  a  net  ion  production/loss  frequency,  vw ,  (which 
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is  independent  of  the  ionization  frequency,  z^),  a  frequency  for  ion  radial  friction,  vr.  and  the 
ion  temperature,  T*m,  at  the  point  separating  inner  and  outer  radial  flows. 

The  plasma  balance  condition  of  Tonks-Langmuir  stated  that  geometrical  magnitudes  de¬ 
termined  Vi  and,  as  a  result,  the  electron  temperature.  In  the  present  model,  the  plasma 
balance  condition  yields  vw  as  a  function  of  the  effective  radial  friction,  the  ionization  rate, 
the  temperature  T*m,  and  geometrical  magnitudes,  whereas  Te  remains  a  free  parameter  (to  be 
determined  from  the  axial  conditions  of  the  discharge).  The  investigation  of  the  plasma  bal¬ 
ance  condition  reveals  two  important  conclusions.  First,  vw.  beyond  being  proportional  to  the 
channel  width,  depends  significantly  on  vr  in  a  way  such  that  a  strong  radial  friction  can  reduce 
significantly  the  loss  of  plasma  current  into  the  lateral  walls.  And  second,  v,t  has  practically  no 
direct  influence  on  vw ;  the  indirect  influence  comes  from  the  contribution  of  i \  to  vr. 

The  equation  for  the  ion  internal  energy  includes  two  source  terms:  one  due  to  the  velocity 
dispersion  of  the  ion  fluid  because  of  ionization,  and  the  second  one  coming  from  the  mean 
temperature  of  the  ion  sources  (ionization  and  axial  flow).  As  a  consequence,  T*m  is  partially 
independent  of  the  temperature  of  the  neutral  gas,  Tn.  Regarding  the  thermodynamic  response 
of  the  ions,  it  is  shown  that:  (i)  an  isothermal  behavior  exists  only  for  the  zero-ionization  and 
zero-temperature  limit;  (ii)  an  isentropic  response  takes  place  for  zero-ionization,  with  radial 
rarefaction  producing  ion  cooling;  and  (iii)  for  Vi  >  0  and  TiM  >  0,  heating  due  to  ionization 
competes  with  the  rarefaction  cooling  in  shaping  the  ion  temperature  profile. 

The  analysis  of  channels  with  different  annular  ratios,  shows  that  asymmetries  between  the 
profiles  of  the  inner  and  outer  presheaths  are  relevant  for  <  0.7,  roughly.  Point  M  of  maximum 
potential  moves  towards  the  inner  wall  as  /3  decreases,  and  plasma  gradients  are  larger  in  the 
inner  presheath.  In  spite  of  that  presheath  asymmetry,  for  v^  —  0  the  potential  drops  and  the 
currents  into  the  walls  are  the  same  in  the  inner  and  outer  presheath.  For  v,t  >  0(1),  smaller 
potential  drops  and  larger  wall  currents  are  found  in  the  inner  presheath.  Plasma  currents 
remain  constant  inside  the  sheaths  and  the  sheath  potential  drops  are  practically  independent 
of  the  presheath  solution. 

To  get  an  idea  of  the  practical  relevance  of  the  different  parameters  controlling  the  presheath 
solution  typical  values  for  a  SPT-100  type  of  Hall  thruster  are  /3  ~  0.7,  T)m  <C  1,  v0  ~ 

2  —  4  x  105s-1,  Vi  ~  1  in  the  ionization  region  and  i>i  <C  1  outside  it,  and  vr  ~  1  in  most  of  the 
channel  [7]. 

Finally,  research  is  pursued  in  two  directions.  First,  ceramic  materials  may  present  a  sec¬ 
ondary  emission  yield  of  order  unity[55,  36],  which  is  not  immediate  to  include  in  a  presheath/sheath 
model  between  two  walls[9].  Second,  to  determine  completely  the  conditions  of  the  plasma  dis¬ 
charge  requires  the  simultaneous  solving  of  the  coupled  radial  and  axial  models[7]. 
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Figure  1.1:  A  generic  radial  section  of  a  long  annular  channel  between  two  dielectric  walls. 
9xi  —  nxivxi  and  gri  —  represent  axial  and  radial  particle  fluxes;  <j)(r )  sketches  the  radial 

profile  of  the  electric  potential,  which  is  maximum  at  point  M.  The  transitions  to  the  sheaths  are 
at  points  Q  and  Q'.  Since  the  presheath/sheath  structure  corresponds  to  the  zero  Debye-length 
limit  one  has  r\y  —  tq  and  r'w  —  r'q. 
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Figure  1.2:  Planar  geometry:  Evolution  of  plasma  parameters  with  vr  for  (T,m,  z'*)  =  (0,  0)  [solid 
line],  (1,  0) [dash-dot  line],  and  (0,  z^r)[dashed  line].  The  dimensionless  particle  flux  is  gr{  — 
(j^iM  VTe/ )  Sri- 
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Figure  1.3:  Planar  geometry:  Net  production  frequency  versus  ionization  frequency  for  = 
0(solid  lines)  and  l(dashed  lines)  and  / vr  =  0.5, 1  and  2  (from  bottom  to  top  for  each  T)m)- 
The  asterisks  correspond  to  the  no-flowing  case,  vw  —  Vi  —  vr. 


Figure  1.4:  Planar  geometry:  Radial  profile  of  ion  temperature  for  =  0.3,  vr  =  Dj.  and  four 
ionization  frequencies:  =  0, 1,  3  and  10. 
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Figure  1.5:  Annular  geometry:  Plasma  profiles  for  different  annular  ratios  (r'Q/rQ  c  0.013, 
0.072,  0.21,  0.43),  T*m  =  0.5,  and  vr  —  Vi  —  0.  Points  Q'  are  at  the  asterisks;  points  M  are  at 
d(f>/ dr | m  —  0  (and  v„  —  0)  separating  the  inner  and  outer  presheaths. 
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Figure  1.6:  Annular  geometry:  Evolution  of  main  plasma  parameters  with  i \  for  =  0.5, 
vr  —  Vi  and  three  annular  ratios,  ft  —  0.1,  0.4  and  0.7.  In  (b)  and  (c)  solid  lines  correspond  to 
the  outer  presheath,  and  dashed  lines  to  the  inner  presheath. 
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Chapter  2 


RADIAL  MODEL  WITH  SECONDARY 
EMISSION 

2.1  INTRODUCTION 

The  transverse  (or  radial)  interaction  of  a  collisionless  plasma  with  two  dielectric  walls  confining 
it  was  treated  originally  by  Tonks  and  Langmuir [54],  They  develop  the  sheath/presheath  theory 
for  the  zero-Debye  length  asymptotic  limit  and  solved  the  fundamental  issues  of  the  problem. 
In  Tonks- Langmuir  model  and  posterior  papers  complementing  it,  [33,  23,  16,  50]  the  ion 
and  electron  losses  to  the  walls  are  balanced  by  volume  ionization.  In  Ch.  1  we  generalized 
that  model  to  plasmas  flowing  longitudinally  between  planar  or  annular  walls.  The  presheath 
problem  depends  then  on  two  new  parameters  related  to  the  axial  flow,  and  it  is  shown  that 
(i)  in  the  absence  (or  excess)  of  ionization,  the  axial  flow  balances  the  radial  losses  of  particles 
to  the  walls,  and  (ii)  axial  gradients  introduce  a  friction  effect  on  the  radial  acceleration  of  the 
ions,  which  leads  to  larger  potential  drops  in  the  presheath  and  lower  losses  to  the  walls. 

All  the  above  models  consider  implicitly  a  passive  (or  non-emitting)  wall.  However,  in  many 
cases  electron  emission  by  the  wall  is  important.  For  certain  ceramic  materials,  electron  impact 
produces  secondary  emission  yields  (the  ratio  between  the  secondary  and  primary  fluxes  at  the 
wall)  above  100%  [43,  36].  Secondary  electrons  interact  with  the  electrostatic  field  and  the 
main  plasma  and  modify  sheath  and  presheath  characteristics.  The  collisionless  sheath  with 
secondary  emission  was  studied  by  Hobbs  and  Wesson [34],  who  demonstrated  that  the  sheath 
becomes  charge-saturated  for  a  secondary  emission  yield  close  to  100%.  However,  they  studied 
the  sheath  exclusively,  and  they  did  not  match  it  to  any  model  for  the  bulk  of  the  plasma. 
To  close  the  sheath  equations,  they  impose  the  usual  quasineutrality  and  sonic  (i.e.  Bohm) 
conditions  at  the  sheath/plasma  boundary. 

Schwager  formulated  a  planar,  kinetic  model  with  secondary  emission  of  the  region  between 
one  emitting  wall  and  a  source-  or  edge-plasma  [51].  That  region  consisted  of  a  (negative)  wall 
sheath,  followed  by  a  quasineutral  region,  and  a  (positive)  sheath  close  to  the  edge-plasma. 
Stephens  and  Ordonez  studied  practically  the  same  model  [52],  and  Jolivet  and  Roussel  apply 
the  last  one  to  the  case  of  a  Hall  thruster  [36].  In  these  three  models,  secondary  electrons  are 
assumed  to  accelerate  freely  across  the  whole  plasma  region  between  the  emission  wall  and  the 
edge-plasma.  The  secondary  population  is  thermalized  at  the  edge-plasma  and  re-emitted  back 
as  primary  electrons  to  the  central  plasma  region.  This  wall  plus  edge-plasma  model  does  not 


23 


reproduce  satisfactorily  the  case  of  a  plasma  confined  between  two  walls,  with  no  intermediate 
edge-plasma.  In  particular,  this  edge-plasma  does  not  act  as  a  ’symmetry  mirror’  in  a  planar 
geometry:  (i)  the  solution  for  the  secondary  population  is  not  symmetric,  since  that  population 
would  be  ’free’  in  half-channel  and  ’thermalized’  in  the  other  half;  and  (ii)  the  edge-plasma 
tends  to  create  a  sheath-like  region  close  to  it  which  is  not  expected  in  the  usual  two- wall  case. 

As  the  works  just  commented  suggest,  the  main  difficulty  lies  in  the  treatment  of  the  sec¬ 
ondary  electrons  in  the  bulk  of  the  plasma.  In  a  recent  conference  we  proposed  that  a  consistent 
presheath/sheath  theory  with  secondary  emission  could  be  developed  by  treating  differently  the 
electron  populations  in  each  region  [9].  The  present  chapter  develops  in  detail  that  model.  The 
two  main  pieces  of  it  will  be  the  Hobbs- Wesson  model  for  the  sheath  (extended  here  to  the  case 
of  warm  ions)  and  our  presheath  model  for  flowing  plasmas  and  two  annular  walls. 

An  application  of  the  resulting  model,  which  is  of  particular  interest  to  us,  is  the  plasma 
discharge  in  a  Hall  thruster  chamber  made  of  ceramic  walls.  The  discharge  is  axisymmetric 
and  particle  and  energy  losses  to  the  walls  are  known  to  affect  greatly  the  thruster  performance 
[46,  26,  38,  7].  Therefore,  a  good  understanding  of  the  radial  plasma- wall  interaction  is  crucial 
to  predict  correctly  the  axial  response  of  the  discharge. 

The  rest  of  the  chapter  is  organized  as  follows.  In  Sec.  II  we  formulate  the  main  hypotheses 
of  the  model.  In  Sec.  Ill  we  discuss  the  sheath  solution.  In  Sec.  IV  we  match  the  presheath 
model  to  the  sheath  solution  and  we  compute  results  for  the  complete  solution.  In  Sec.  V  we 
discuss  the  main  results. 


2.2  MODEL  FORMULATION 

Figure  2.1(a)  sketches  the  radial  model.  The  plasma  flows  between  two  annular  walls  placed, 
in  cylindrical  coordinates  ( x,r,9 ),  at  r  —  rw1  and  r  =  rw  =  +  hc,  with  hc  the  channel 

width.  Here,  we  are  interested  in  the  radial  structure  of  the  plasma  discharge  at  a  given  axial 
position  x.  For  dielectric  walls,  the  zero  current  condition  implies,  in  general,  the  potential 
profile  cf)(r )  to  present  a  maximum  at  a  central  point  of  the  channel  [point  M  in  Fig.  2.1(a)] 
and  to  decrease  monotonically  towards  the  two  walls  [points  W  and  W’].  Except  for  very  large 
voltages,  secondary  electron  emission  is  due  mainly  to  the  impact  of  those  electrons  from  the 
main  population  with  enough  energy  to  reach  the  walls.  Depending  on  their  posterior  evolution 
inside  the  plasma,  secondary  electrons  can  be  divided  in  two  basic  types:  (i)  a  ’free’  group 
constituted  by  those  crossing  the  plasma  from  the  emission  wall  and  being  collected  back  at 
the  opposite  wall,  without  interacting  with  other  particles,  and  (ii)  a  ’trapped’  group  formed 
by  those  suffering  occasional  collisions  with  other  particles  and  multiple  electrostatic  reflexions 
within  the  confining  potential  </>(r),  until  they  thermalize  with  the  bulk  electron  population. 

The  ’free’  secondary  population  is  simple  to  model;  on  the  contrary,  at  each  instant  of  time, 
there  are  ’trapped’  secondary  electrons  with  different  temporal  histories.  The  crucial  issue  of 
any  stationary  model  is  to  figure  out  which  is  the  temporally-averaged  collective  distribution 
of  this  trapped  population.  Our  proposal  is  to  take  advantage  of  the  presheath/sheath  theory 
and  of  a  macroscopic  formulation  for  a  dual  treatment  of  the  ’trapped’  secondary  electrons. 
The  presheath/sheath  theory  applies  to  the  zero  Debye-length  limit, 
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with  \d  the  Debye  length,  and  A coi  the  relevant  mean  free-path  for  collisional  processes.  The 
advantage  of  using  this  two-scale  analysis  is  based  on  three  features:  (i)  the  sheaths  are  treated 
as  surface  discontinuities  in  the  quasineutral  scale,  but  as  as  collisionless  and  ionization-free 
regions  in  their  inner  scale;  (ii)  presheath  and  sheaths  are  solved  separately  in  their  own  scales, 
and  then  matched  at  the  presheath/sheath  transition;  and  (iii)  the  matching  or  transition 
condition  is  unique  and,  for  a  plasma  with  zero  radial  flow  somewhere  in  the  presheath,  consists 
of  the  sonic  Bohm  condition  [40,  3]. 

The  model  we  propose  for  the  plasma  is  based  in  the  following  hypotheses: 

(i)  Within  the  thin,  collisionless  sheaths  there  are  two  independent  electron  populations: 
population  p  of  primary,  confined  electrons  coming  from  the  bulk  of  the  plasma  and  population 
s  of  secondary,  wall-emitted  electrons.  Population  s  is  accelerated  freely  by  the  sheath  potential 
from  the  wall  to  the  presheath  boundary. 

(ii)  We  consider  that  A coi  <C  hc,  so  that  all  secondary  electrons  remain  trapped  in  the  bulk 
of  the  plasma.  [This  restricts  the  present  model  to  the  total-trapping  case.] 

(iii)  Within  the  presheath,  all  (trapped)  electrons,  irrespective  of  their  origin  (ionization, 
axial  flow,  or  wall-emission),  are  grouped  in  a  single  population  e. 

(iv)  There  is  a  unique  population  (i)  of  ions  in  sheath  and  presheath.  Quasineutrality  in 
the  presheath  means  rq  ~  ne  at  every  point  (and  not  just  globally,  as  Ref.  [52]  states). 

(v)  To  match  consistently  the  solutions  of  presheath  and  sheaths  three  basic  laws  are  im¬ 
posed  at  each  presheath/sheath  transition:  the  conservation  of  the  total  electron  density,  the 
conservation  of  the  net  electron  flux,  and  the  continuity  of  the  ion  flux. 


2.3  SHEATH  ANALYSIS 

2.3.1  Sheath  solution 

In  general,  we  will  refer  the  discussion  only  to  the  sheath  tied  to  the  outer  wall;  the  inner 
sheath  behaves  in  the  same  way.  The  outer  sheath  extends  from  the  the  wall  (point  W)  to  the 
presheath  transition  (point  Q).  If  the  sheath  is  not  charge-saturated  the  electric  potential  is 
monotonic  with  (f>wQ  —  ~  <fov  >  0-  In  the  zero  Debye-length  limit  of  Eq.  (2.1),  the  sheath 

satisfies  the  conservation  equations, 


riiVri  —  const  =  gri: 

(2.2) 

^171^  +  +  ecf)  =  const, 

(2.3) 

T,; /nf'A  —  const, 

(2.4) 

ecf)  —  e(j)w 

np  =  npW  exp 

J-p 

(2.5) 

TisVrs  =  const  =  grs, 

(2.6) 

-mev^s  —  ecf)  —  const, 
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the  flux  conditions, 


9ri  9rp  T  9rsi 

(2.8) 

9rp  —  TlpW  ^Tp /27T777.e, 

(2,9) 

9rs  (Tp  )  9rp  ■ 

(2,10) 

and  the  (quasiplanar)  Poisson  equation, 

eo  d2(f) 

,9  =np  +  ns  rii- 
e  drz 

(2,11) 

The  constants  for  the  ion  equations  (2.2)-(2.4)  are  set  at  point  Q  and  must  match  the  values 
obtained  at  the  presheath  side.  Since  the  wall  temperature,  Tw.  is  generally  much  smaller  than 
the  plasma  temperature  Tp,  the  mean  energy  of  the  secondary  electrons  when  they  leave  the 
wall  can  be  neglected  (except  in  the  very  vicinity  of  point  W)  in  Eq.  (2.7).  Equation  (2.8)  is  the 
zero-current  condition,  and  Eq.  (2.9)  corresponds  to  assume  a  quasi-Maxwellian  distribution 
for  population  p. 

In  Eq.  (2.10),  the  effective  secondary  emission  yield  of  the  wall,  SW(TP),  depends  on  both 
Tp  and  the  wall  material.  For  the  present  purposes,  a  good-enough  expression  of  the  emission 
yield  produced  by  a  monoenergetic  beam  of  energy  E  per  electron  is  dw(E)  ~  y/£’/£'° ,  with 
E °  the  temperature  for  100%  yield  (in  the  range  E°  ~  10  —  100V,  say) [26,  36].  Averaging  over 
the  quasi-Maxwellian  population  p  the  effective  yield  satisfies 

Sw(Tp)  -  \Jtp/Ew,  (2.12) 

with  Ew  ~  0.57 E^.  This  last  relationship  allows  us  to  use  Ew/Tp  instead  of  Sw  as  input 
parameter  for  the  sheath  problem. 

Equation  (2.11)  states  that 


c  = 


r  ~  rw 
A  d 


with  A  h  — 


eo  Tp 
e2niQ ' 


(2.13) 


is  the  natural  spatial  variable  in  the  sheath.  The  integration  of  Poisson  equation  (2.11)  together 
with  Eqs.(2.4)-(2.7)  yields  the  conservation  equation 


d(f)V2  fd(j)  y 
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with 


U(e<f>)  =  ripTp  +  nsmev2s  +  nt(rntv2t  +  T{) 
the  Sagdeev’s  potential [45].  First  and  second  derivatives  of  this  function  are 

U'{e<f> )  =  np  +  ns  -  nt, 
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Equation  (2.14)  states  the  balance  of  the  dynamic  and  electric  pressures  of  the  whole  plasma 
across  the  sheath.  There  are  three  boundary  conditions  related  to  Poisson  equation.  First, 
plasma  quasineutrality  at  point  Q  means 


(2.17) 


Second,  for  a  solution  to  exist  around  point  Q  with  \d(j)/d(\  increasing  towards  the  wall, 
Eq.  (2.14)  requires  that  Uq  >  0.  In  terms  of  the  ion  velocity  (and  neglecting  the  case 
mivriQ  <  fTjQ  which  does  not  lead  to  a  valid  solution  at  the  wall)  this  condition  is  equiv¬ 
alent  to 


2  ^  5 

miVriQ  >  -TiC 
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n'pQ  +  nSQ 


Tl pQ  W'sQ 


(2.18) 


which  is  the  sonic/supersonic  Bohm  condition  for  this  three-species  sheath.  Since  the  quasineu¬ 
tral  solution  of  the  presheath  accepts  only  a  subsonic/sonic  transition  (see  Sec.  2.4.2),  the  sonic 
condition 

Uq  =  0,  (2-19) 

is  the  correct  one  for  the  sheath/presheath  transition.  Finally,  the  analysis  of  the  transition 
region  around  point  Q  shows  that  the  electric  field,  when  measured  in  the  sheath  scale,  tends 
asymptotically  to  zero  at  the  sheath  side  of  point  Q,[40,  48] 
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The  preceding  equations  and  conditions  show  that  the  dimensionless  solution  of  the  sheath 
depends  on  just  two  parameters,  5W  (or  Ew/Tp )  and  T }q  —  TiQ/Tp.  Hobbs  and  Wesson  treated 
the  case  of  cold  ions,  T^q  =  0.  Figure  2(a)  shows  potential  profiles  for  different  emission  yields 
and  two  ion  temperatures.  For  a  given  potential,  a  warm  ion  population  has  a  larger  density 
and  a  lower  velocity  than  a  cold  one,  according  to  Eq.  (2.3);  this  implies  a  smaller  Debye  length 


2.3.2  The  charge-saturated  regime(CSR) 


Hobbs  and  Wesson  showed  (for  %q  —  0)  that  the  sheath  reaches  the  charge  saturation  limit  (CSL) 


at  the  wall,  characterized  by  d(ecf)/Tp)/d( 
drop  are 


w 


Sw  =  5*w  ~  1  -  8.3 


0,  when  the  secondary  yield  and  the  potential 

e(f*WQ 


1 t>wQ  —  1-02. 


(2.22) 


For  warm  ions  (i.e.  T*q  >  0),  the  situation  is  similar.  We  find  that  the  ion  temperature  has  an 
irrelevant  influence  on  the  value  of  Sw  at  the  CSL:  for  xenon,  one  has  6*w  ~  0.983  for  Tjq  =  0 
and  5 ^  ~  0.981  for  T \q  =  1.  The  effect  of  Tjq  on  other  magnitudes  is  small  but  not  negligible: 
the  sheath  potential  drop  at  the  CSL,  </>^Q(TjQ),  decreases  about  a  17%  when  Tjg  goes  from  0 
to  1. 

For  5W  >  5 ^  the  sheath  is  in  the  charge-saturated  regime(CSR).  A  potential  well  is  formed 
near  the  wall  to  turn  back  the  excess  of  secondary  flux.  Its  magnitude,  A <f>w,  must  be  such  that 
the  secondary  flux  satisfies  the  CSL  condition  grs  =  —d^g^  at  the  point  of  minimum  potential. 
A  simple  estimate  yields 

A(f>w  ~  (Tw/e)\n(5w/5*w).  (2.23) 


Then,  in  the  asymptotic  limit  Tw/Tp  — >  0,  the  potential  well  can  be  neglected  and  the  dimen¬ 
sionless  solution  for  the  whole  charge-saturated  regime  (CSR)  coincides  with  that  of  the  CSL. 
In  particular,  one  has,  for  the  CSR, 


=  const  =  4>*WQ(fiQ). 


(2.24) 


A  final  observation  is  that  Eq.(2.22)  implies  that  gri/\grs\  —  8.3 ^/me/mj,  a  result  which 
agrees  in  order  of  magnitude  with  the  well-known  Langmuir  condition,  gri/\grs\  —  me/mi1 
for  a  classical  strong  double  layer  (characterized  by  e(f>wQ/Tp  1)  [40]. 


2.4  COMPLETE  SOLUTION 

2.4.1  Presheath  solution 

The  solution  of  the  sheath  is  universal  for  any  value  of  the  plasma  density  and  the  electron 
temperature  at  the  sheath  entrance.  In  particular,  the  net  plasma  current  to  the  wall,  which  is 
constant  across  the  sheath,  is  determined  by  the  presheath. 

For  total  trapping  of  secondary  emission,  the  quasineutral  plasma  of  the  presheath  consists 
of  (i)  the  population  i  of  ions,  which  accelerate  from  point  M  of  maximum  potential  and  zero  ion 
velocity  towards  the  sheath  transition  points  Q’  and  Q,  and  (ii)  a  single  population  e  of  confined 
electrons,  which  are  in  Maxwell- Boltzmann  equilibrium  at  temperature  Te.  Since  population  e 
differs  from  population  p  of  the  sheath,  their  temperatures  (which  express  the  average  random 
energy  of  the  two  populations)  are  different,  in  general,  and  the  relation  between  them  is  part 
of  the  complete  solution. 

According  to  Ch.  1,  for  the  general  case  of  a  plasma  flowing  axially  between  two  annular 
walls,  the  presheath  problem  depends  on  the  channel  geometry  (annular  ratio  /3  —  rw'/rw ,  and 
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channel  width  hc) ,  the  plasma  conditions  (ion  density  and  temperature  at  point  M,  njM  and 
TiM,  and  electron  temperature  Te),  and  effective  frequencies  for  net  plasma  production  (vw), 
ionization  (v,),  and  radial  friction  (vr).  The  presheath  problem  states  that  the  transition  to 
the  (outer)  sheath  is  defined  by  the  sonic  Bohm  condition  [40,  3] 

mivriQ  =  ^ TiQ  +  Te.  (2.25) 

The  dimensionless  presheath  solution  depends  on  just  four  dimensionless  parameters:  B.  TiM  = 
TiM/Tei  Di  =  Vi/v o,  and  vr  =  vr/v 0,  with  =  h~x  ^Te/nii.  The  relationship 

vw  =  vw{P,  TiM,  Di,  vr)  (2.26) 

is  known  as  the  plasma  balance  equation;  for  the  no-flowing  case  of  Tonks-Langmuir,  one  has 
vw  —  Vi  —  vr  and  the  plasma  balance  equation  becomes  Vi(B-,  Tm)-  Figure  2.5,  commented 
below,  illustrates  the  presheath  solution  with  a  particular  example. 

2.4.2  Presheath/sheath  matching 

The  correct  matching  of  the  presheath  and  sheath  solutions  comes  out  from  three  basic  con¬ 
servation  equations,  which,  indeed,  have  been  imposed  already:  first,  there  is  the  particle 
conservation, 

npQ  +  n sQ  =  ^eQ  =  niQ'-, 

then,  there  is  the  continuity  of  the  electron  flux, 


9rpQ  T  9rsQ  —  9reQ  —  9r%Qi 


Fig.  2.1(b);  and,  third,  there  is  the  continuity  and  compatibility  of  the  ion  flux,  which  requires 
to  satisfy  the  sonic  Bohm  conditions  at  the  sheath  and  presheath  sides,  Eqs.  (2.19)  and  (2.25). 

Of  course,  for  no  secondary  emission,  populations  p  and  e  coincide  and  Eqs.  (2.19)  and 
(2.25)  are  the  same  condition.  The  demonstration  that  this  sonic  condition  yields  the  only  valid 
transition  is  well  established  [54,  48,  3].  For  5W  >  0,  to  consider  that  Eqs.  (2.19)  and  (2.25) 
are  the  same  condition  is  (partially,  at  least)  a  postulate  based  in  fundamental  properties  of 
the  ion  flow  for  the  well-known  zero-emission  case.  [A  complete  demonstration  of  the  postulate 
would  require  to  specify  the  transformation  of  populations  p  and  s  into  population  e  in  a 
transition  layer  centered  at  point  Q.]  Adding  support  to  it,  this  postulate  or  condition  closes 
the  presheath/sheath  matching  since  it  provides  the  relationship  we  needed  between  Tp  and  Te: 
equating  Eqs.  (2.19)  and  (2.25)  one  has 

^ _ TpnSQ 

Tp  _  2e(f)wQnpQ  ^  1 

Te  “  !  +  ^ 

npQ 


This  way  to  relate  the  two  temperatures  is  very  adequate  to  our  macroscopic  model,  more  than 
a  method  based  on  the  electron  distribution  functions  at  the  two  sides  of  point  Q.  Figure  2.4 
shows  the  evolution  of  the  temperature  ratio  Tp/Te  with  5W.  which  is  very  similar  to  that  of 
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ripQ/niQ  in  Fig.  2.3(c).  The  fact  that  both  are  close  to  one  even  at  the  CSR,  leads  us  to  conclude 
that  the  details  of  the  modeling  of  the  electron  populations  is  not  very  critical  for  the  overall 
model. 

In  a  recent  two-dimensional  model  of  the  plasma  discharge  in  a  Hall  thruster,  Keidar  et  al.  use 
a  radial  presheath/sheath  model  with  secondary  emission  [37].  Their  formulation  is  similar  to 
ours  for  cold-ions  (T)  =  0)  but  the  presheath/sheath  matching  is  solved  unsatisfactorily.  First, 
they  do  not  explain  how  they  relate  the  trapped  populations  in  sheath  and  presheaths  and  just 
take  Tp  —  Te,  even  at  the  CSR.  Second  and  more  important,  they  impose  a  subsonic  ion  flow  at 
the  sheath  entrance,  violating  the  supersonic/sonic  condition  (2.18).  The  authors  do  not  show 
radial  profiles  to  support  their  hypotheses. 

2.4.3  Results 

Figure  2.5  shows  an  example  of  a  matched  sheath/presheath  solution  for  an  annular  channel. 
Profiles  are  shown  in  the  presheath  scale,  where  the  sheaths  are  two  wall  discontinuities;  on  the 
contrary,  the  sheaths  are  semi-infinite  in  the  inner  scale  ((q  =  — oo  in  Fig.  2.2).  Wall  conditions 
in  Fig.  2.5  are  such  that  the  two  sheaths  are  in  the  CSR.  The  profiles  of  the  inner  and  outer 
presheaths  are  asymmetrical:  larger  potential  drops  and  lower  currents  are  found  in  the  outer 
presheath.  The  presheath  asymmetry  introduces  a  very  small  asymmetry  between  the  two 
sheaths  (less  than  1%).  An  interesting  observation  is  that  the  potential  drop  in  the  presheath 
is  not  negligible  at  all;  indeed,  in  Fig.  2.5,  the  potential  drop  in  the  outer  presheath  is  larger 
than  in  the  sheaths,  4>qm  >  4>wq ■  The  profile  of  the  ion  temperature  illustrates  the  competition 
between  the  heating  effect  due  to  ionization  and  the  cooling  effect  due  to  rarefaction,  this  last 
one  dominating  totally  in  the  sheaths  and  justifying  the  small  effect  of  the  ion  temperature  on 
the  sheath  solution. 

The  sheath  solution  showed  that  the  relative  potential  drop  ecfrwQ/Tp  decreases  as  Sw  in¬ 
creases,  and  remains  constant  in  the  CSR.  Figure  2.6(a)  shows,  for  different  wall  materials,  the 
evolution  of  the  sheath  potential  drop  with  the  electron  temperature  in  the  bulk  of  the  channel, 
Te.  There  is  a  relative  maximum  of  (f>wq{Te ),  followed  by  a  minimum  just  at  the  CSL,  and  then 
a  linear  increase  with  Te:  one  has  e^q/Tg  ~  0.88  (for  xenon  and  =  0).  An  important 
point  is  that  4>wq  is  independent  of  the  wall  material  once  in  the  CSR.  The  potential  drops  in 
the  presheaths,  cj)qM  and  4>q'Mi  must  be  computed  from  the  presheath  solution  of  Ch.  1. 

At  a  given  radial  cross-section,  the  energy  fluxes  of  ion  and  electron  into  the  (outer)  wall 
are 

QiW  —  9riQ  (v^iW  +  Vxi)  +  §Tjiy)  =  driQ^^WQ  +  \Te  +  Y^iQ  +  i 

(2.28) 

where  vxi  is  the  (average)  velocity  of  the  ion  axial  flow  in  that  radial  section.  The  first  observed 
feature  is  that  electron  losses  depend  strongly  on  the  secondary-emission.  Indeed,  if  there 
is  no  axial  flow,  ion  energy  losses  are  about  3  times  larger  than  electron  energy  losses  for 
5W  —  0,  whereas  they  become  equal  for  5W  ~  60  —  70%,  and  electron  energy  losses  become  two 
orders  of  magnitude  larger  than  ion  ones  in  the  CSR.  Warm  and  axial  flow  effects  increase  ion 
energy  losses:  in  the  acceleration  region  of  a  SPT-100  Hall  thruster  and  for  charge-saturated 


QeW  —  drpQ^Tp  —  QriQ  ^ 


-S' 
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conditions,  this  increase  can  be  up  to  one  order  of  magnitude.  Figures  6(b)-6(c)  show  the 
evolution  of  electron  and  ion  energy  losses  with  Ew  and  Te.  Notice  the  large  gradients  of  Qew 
with  Te  below  the  SCL.  Again,  the  influence  of  the  wall  material  disappears  at  the  CSR,  when 
Qew  scales  as  T^2 . 


2.5  DISCUSSION 

The  presheath/sheath  model  with  total  trapping  of  secondary  emission  presented  here  is  based 
in  considering  wall-emitted  electrons  as  a  free  population  within  the  thin,  collisionless  sheaths, 
and  as  part  of  a  unique  population  of  confined  electrons  within  the  presheath.  As  a  consequence 
the  populations  of  confined  electrons  in  sheath  and  presheath  are  somewhat  different.  The 
presheath  and  the  sheaths  are  solved  separately  and  then  matched  using  continuity  equations  for 
the  plasmas  density  and  fluxes.  The  sheath  model  is  an  extension  of  the  Hobbs- Wesson  model  to 
warm  ions,  and  presents  a  charge-saturated  regime.  The  presheath  model  is  the  one  developed 
recently  by  us  for  plasmas  flowing  axially,  and  of  application  to  planar,  cylindrical,  and  annular 
geometries.  The  imposition  of  sonic  Bohm  conditions  at  the  two  sides  of  the  transition  point 
closes  the  problem  and  relates  the  temperature  of  the  confined  electron  populations  in  sheaths 
and  presheath.  It  is  found  that,  even  at  charge  saturation,  the  density  of  secondary  electrons 
entering  into  the  presheath  is  small  and  the  temperatures  of  the  confined  electron  populations 
in  sheaths  and  presheath  are  similar. 

The  complete  sheath/presheath  solution  depends  on  five  dimensionless  parameters:  0. 
ur,  Tm,  and  Ew/Te.  The  first  four  determine  the  presheath  solution  whereas  the  sheath  so¬ 
lution  depends  mainly  on  the  fifth  one.  At  the  charge-saturated  regime  the  solution  becomes 
independent  of  the  wall  material.  A  relevant  feature  is  that  particle  losses  to  the  walls  are 
determined  totally  by  the  presheath,  whereas  the  sheath  solution  gives  only  the  energy  losses 
per  particle.  Electron  energy  losses  change  drastically  with  secondary  emission.  Ion  energy 
losses  depend  weakly  on  Sw.  have  a  contribution  from  the  plasma  axial  flow,  and  dominate  for 
low  secondary  emission. 

The  sheath  does  not  vanish  at  the  charge-saturated  regime.  The  sheath  potential  drop 
presents  a  maximum  for  an  intermediate  value  of  5W  and  a  relative  minimum  at  the  charge- 
saturated  limit;  then,  it  becomes  proportional  to  the  electron  temperature.  A  non-zero  ion 
temperature  does  not  modify  the  emission  yield  necessary  for  charge  saturation  but  diminishes 
the  relative  sheath  potential  at  the  charge-saturation  regime.  The  potential  drops  in  the  inner 
and  outer  presheaths  are  negligible  never.  Furthermore,  they  can  be  larger  than  the  sheath 
drops  for  large  secondary  emissions  and  large  radial  frictions.  For  annular  channels,  the  inner 
and  outer  presheaths  are  asymmetric,  but  the  two  sheaths  are  practically  identical  (for  the 
same  wall  material). 

The  model  developed  here  is  for  total  trapping  of  the  secondary  emission,  which  corresponds 
to  the  collisional  regime,  A^  <C  A coi  <C  hc.  For  the  collisionless  regime,  A coi  hc.  which  is  more 
typical  of  Hall  thruster  discharges,  secondary  electrons  could  cross  the  whole  presheath  from 
the  emission  wall  and  be  collected  at  the  opposite  wall  if  they  had  enough  radial  energy  to  reach 
it.  For  instance,  for  the  potential  profile  of  Fig.  2.5(a),  A coi  hc,  and  cold  walls,  electrons 
emitted  from  the  outer  wall  will  be  collected  back  at  the  inner  wall,  whereas  electrons  emitted 
from  the  inner  wall  will  be  reflected  back  at  the  outer  wall  and  will  become  either  trapped  in 
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the  plasma  or  collected  back  by  the  inner  wall  too.  However,  the  presence  of  free  populations  of 
secondary  electrons  modifies  in  turn  the  radial  of  the  electric  potential  in  sheaths  and  presheath. 
Therefore, the  present  model  is  not  applicable  to  partial  trapping  of  secondary  emission  and  a 
more  general  must  be  developed. 
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Figure  2.1:  (a)  A  generic  radial  section  of  a  long  annular  channel  between  two  dielectric  walls; 
gx  and  gr  represent  axial  and  radial  particle  fluxes,  respectively;  <j)(r )  sketches  the  electric  po¬ 
tential,  which  is  maximum  at  at  intermediate  point  M;  points  Q  and  Q'  represent  the  transitions 
to  the  two  sheaths,  which  actually  are  discontinuities  in  the  quasineutral  scale,  (b)  Treatment 
of  the  electron  populations  in  radial  presheath  and  sheath:  in  each  sheath,  primary(p)  and 
secondary (s)  electron  populations  are  defined;  in  the  presheath,  a  single  population (e)  is  de¬ 
fined;  zero-current  condition  plus  conservation  conditions  at  the  presheath/sheath  transition 
mean  9ri  —  SVe  —  9rp  T  grs  there. 
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Figure  2.2:  (a)  Potential  profile  for  different  values  of  the  secondary  emission  yield  and  TiQ  =0 
(solid  lines)  and  1  (dashed  lines).  The  cases  Sw  =  5*m (Ttq )  correspond  to  d<p/d(\w  =  0,  marking 
the  charge-saturation  limit  of  the  sheath,  (b)  Profiles  of  the  different  plasma  species  for  T%q  =0 
and  5W  =  5^;  fii  —  rip  —  ns  is  the  net  electric  charge  in  the  sheath. 
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Figure  2.3:  Evolution  with  5W  of  the  sheath  potential  drop,  the  electric  held  at  the  wall,  and  the 
density  of  secondary  electrons  at  point  Q  for  T)q  =0  (solid  lines)  and  1  (dashed  lines).  Circles 
correspond  to  the  charge  saturation  limit.  At  the  SCL  <P*Wq  goes  from  1.02  at  Tjg  =  0  to  0.846 
at  TiQ  =  1. 
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Figure  2.4:  Evolution  with  8W  of  the  temperature  ratio  for  the  confined  electron  populations  p 
of  sheath  and  e  of  presheath  for  %q  =0  (solid  lines)  and  1  (dashed  lines).  At  the  SCL  (Tp/Te)* 
goes  from  0.865  at  =  0  to  0.819  at  Ttq  =  1. 
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Figure  2.5:  Complete  presheath/sheaths  solutions  for  an  annular  channel  with  /3  =  0.285, 
TiM/Te  =  0.50,  vr  =  Vi  =  3.49,  and  charge-saturated  conditions  at  the  two  sheaths,  Q'W'  and 
QW .  The  representation  is  in  the  quasineutral  scale.  In  (c)  g„  =  (riiM \/Te/ rrii)-1  gri  is  the 
dimensionless  flux,  which  is  constant  across  each  sheath. 
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Figure  2.6:  (a)  Sheath  potential  drop,  (b)  electron  energy  losses,  and  (c)  ion  energy  losses  as 
function  of  Te  for  xenon,  Tiq/Te  =  0,  and  three  wall  materials:  Ew  =  15,  25  and  40  eV.  The 
curves  are  independent  of  Ew  at  the  CSR.  Plots  (b)  and  (c)  are  for  =  1018m-3;  for  other 
densities,  particle  and  energy  losses  scale  proportional  to  ntQ. 
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Chapter  3 

PAixTIAL  TixAPPING  OF 
SECONDAixY  ELECTixON  EMISSION 


3.1  INTRODUCTION 

Chapter  2  presented  a  macroscopic  model  of  the  interaction  of  a  plasma  with  two  parallel, 
electron-emitting,  dielectric  walls.  The  model  is  of  application  to  quiescent  and  flowing  plasmas, 
and  planar,  annular,  and  cylindrical  geometries.  The  zero  Debye-length  limit  was  applied  to 
develop  a  sheath/presheath/sheath  model  in  the  radial  direction  (i.e.  perpendicular  to  the 
walls),  with  the  maximum  of  the  electrostatic  potential  around  the  center  of  the  channel  cavity. 
From  each  wall,  secondary  electron  emission  (SEE)  is  accelerated  by  the  adjacent  sheath  into 
the  bulk  of  the  plasma.  The  central  issue  of  that  chapter  was  the  treatment  of  the  SEE 
once  in  the  presheath.  The  analysis  of  the  chapter  was  restricted  to  the  collisional  (or  total- 
thermalization  or  total-trapping)  case,  when  the  two  SEE  beams  collide  with  the  main  plasma, 
their  radial  momentum  is  dispersed  in  all  directions,  they  cannot  reach  the  opposite  wall,  and 
they  thermalize  eventually  with  the  main  population  of  confined  electrons.  The  model  took 
advantage  of  the  two-scale  analysis  to  treat  SEE  as  independent  beams  within  the  emission 
sheaths,  but  as  parts  of  a  single  electron  population  within  the  presheath.  The  warm-ion 
extension  of  the  three-species  model  of  Hobbs  and  Wesson  [34]  was  used  in  the  sheaths,  and  a 
two-species  model  of  the  presheath  for  flowing  plasmas  in  annular  geometries  developed  recently 
too  was  used  in  the  presheath  [3].  Appropriate  conditions  at  the  presheath/sheath  transitions 
had  to  be  defined  to  match  consistently  the  different  plasma  models  in  sheaths  and  presheath. 

This  chapter  addresses  the  extension  of  the  model  of  Ch.  2  to  partially-collisional  and 
collisionless  cases.  This  study  is  of  direct  application  to  Hall  thrusters  discharges  [36].  For 
instance  for  a  conventional,  1.5kW  thruster  (like  a  SPT-100)  the  channel  width  is  about  15 
mm,  whereas  the  effective  mean-free  path  for  the  SEE  is  estimated  to  range  between  20  and 
200mm,  depending  on  the  chamber  zone  [49]. 

There  are  several  conclusions  of  the  previous  chapter  which  are  worth  to  remind  now.  First, 
the  losses  of  plasma  current  to  the  walls  are  determined  totally  by  the  presheath  solution, 
whereas  the  energy  lost  by  impacting  electron/ion  depends  on  the  SEE  properties  of  the  walls. 
Second,  for  total-trapping  the  sheaths  become  charge  saturated  for  a  SEE  yield  near  100%. 
Third,  the  presheath  structure  depends  on  several  parameters  and  can  present  larger  potential 
drops  than  the  sheaths.  Four,  for  annular  chambers,  the  presheath  solution  is  asymmetric,  and 
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the  two  sheaths  are  different,  in  general. 

The  extension  to  the  partial-trapping  case  we  present  here  is  however  restricted  to  planar 
chambers  when  the  solution  is  symmetrical  with  respect  to  the  chamber  mid-plane.  The  reason 
is  that  for  annular  chambers,  the  number  of  variables  and  types  of  solutions  increases,  the 
effects  of  partial-thermalization  become  mixed  with  those  due  to  asymmetries,  making  difficult 
a  complete  and  clear  discussion  of  the  problem. 

3.2  MODEL  FORMULATION 

A  sketch  of  the  model  is  shown  in  Fig.  3.1.  The  hypotheses  configuring  it  are  the  following: 

(1)  Zero  net-current  is  considered  at  the  lateral  walls.  This  implies,  in  general,  a  radial 
potential  profile  with  a  single  maximum  at  a  point  M,  placed  at  the  center  of  the  channel  (for 
planar  geometry).  Small  potential  wells  near  the  wall  at  charge-saturation  and  other  possible 
profiles  which  could  arise  for  large  SEE  will  be  commented  later. 

(ii)  We  consider  a  planar  chamber  with  the  the  two  walls  made  of  the  same  material.  Then, 
symmetry  considerations  are  of  application  to  the  plasma  populations  and  the  electric  potential. 

(iii)  The  zero-Debye-length  asymptotic  limit 

^d/hc,  ^d/ Koi 0  (3.1) 

holds,  where  A^  is  the  plasma  typical  Debye  length  hc  the  channel  width,  and  A coi  the  effec¬ 
tive  mean-free-path  for  SEE.  This  allows  us  to  carry  out  a  two-scale  sheath/presheath/sheath 
analysis. 

(iv)  The  plasma  is  constituted  by:  ions(f)  being  attracted  by  the  low  wall  potential;  confined, 
primary  electrons(p);  and  two  SEE  populations(s),  one  from  each  wall.  SEE  from  one  wall  is 
accelerated  away  from  it  by  the  electrostatic  potential  of  the  adjacent  sheath.  Assuming  the 
wall  temperature,  Tw,  much  smaller  than  the  sheath  potential  (which  is  typically  of  the  order 
of  the  electron  temperature  Tp )  the  SEE  enters  the  presheath  as  a  monoenergetic  beam. 

(v)  The  evolution  of  the  SEE  beams  in  the  presheath  depends  on  the  collision  probability, 

5C  =  exp(— A  coi/hc).  (3.2) 

The  collisional  (total-trapping)  case  studied  in  Ch.  2  corresponds  to  5C  ~  1.  In  the  collisionless 
(or  zero-trapping)  limit,  dc  ~  0,  the  SEE  beam  emitted  from  one  wall  is  expected  to  cross  freely 
the  presheath  and  to  be  collected  at  the  opposite  wall,  unless  the  potential  difference  between 
the  two  walls  turn  the  beam  back;  in  this  last  case,  it  will  be  collected  back  by  the  emission  wall. 
Finally,  for  the  partially-collisional  case,  it  is  expected  that  one  fraction  of  each  SEE  beam  be 
trapped  and  thermalized  within  the  main  plasma,  and  the  other  part  be  collected  back  at  the 
wall  with  the  largest  potential. 

(vi)  The  collision  probability  5C  is  used  to  split  each  SEE  beam  in  two  parts:  one,  propor¬ 
tional  to  5C,  which  becomes  trapped  and  thermalized  with  the  main  electron  population,  and 
another,  proportional  to  1  —  5C,  which  crosses  freely  the  presheath,  arrives  to  the  second  wall 
with  exactly  the  (small)  radial  momentum  it  had  at  the  emission  wall,  and  is  collected  by  it.  As 
in  Ch.  2,  in  the  presheath,  the  trapped  SEE  is  included  into  the  confined  electron  population. 
The  free  SEE  is  modeled  in  the  same  way  in  presheath  and  sheath  and,  due  to  its  small  impact 
energy,  does  not  produce  extra  SEE  when  it  is  collected  back. 
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(vii)  The  ion  population  is  considered  cold:  T*/Te  — *  0.  This  is  just  a  minor  simplification 
and  not  a  real  restriction,  since  warm-ion  effects  are  qualitatively  the  same  independently  of 
the  level  of  SEE  trapping.  By  eliminating  the  energy  equation  for  ions  we  focus  the  discussion 
on  the  main  issue  of  the  chapter:  the  plasma  response  with  SEE  partial-trapping. 

3.3  SHEATH  ANALYSIS 

Figure  3.1(b)  sketches  the  different  plasma  populations  in  one  sheath:  primary  electrons(p) 
diffusing  from  the  presheath,  secondary  electrons(s)  emitted  by  the  wall,  and  the  free  secondary 
electrons(/)  remaining  from  the  beam  emitted  by  the  opposite  wall.  The  symmetry  of  the 
problem  means  that  the  conditions  of  populations  /  and  s  are  the  same  except  for  the  direction 
of  their  radial  velocity. 

The  plasma  conservation  equations  for  the  collisionless  sheath  are, 


navra  =  gra  =  const,  (a  =  i,p,  s,  /), 

(3.3) 

e(<t>  —  4>w ) 

np  =  npW  exp 

1p 

(3,4) 

Vri  =  \jv2riQ  + 

(3.5) 

Vrf  =  ~Vt.  =  \jv2rsW  +  ~  M, 

(3,6) 

In  the  last  equation  we  will  take  mev2sW/Te  ~  Tw/Te  ~  0. 

The  ion  flux  into  the  sheath  griq  is  determined  mainly  by  the  presheath.  Assuming  popu¬ 
lation  p  to  be  Maxwellian-like  at  temperature  Tp.  the  flux  of  primary  electrons  at  the  wall  is 
given  by 


9rpQ  WpW  Y  Tp/27T7Tlei  (3.7) 

The  equations  relating  the  different  fluxes  are 

9riQ  =  9rpQ  +  9rsW  +  9rfQi  (3.8) 

9rsW  —  $ w9rpQ j  (3-9) 

9rfQ  =  (1  —  <5C)(  —  9rsw)i  (3.10) 


with  5w(Tp )  the  SEE  yield,  which  depends  on  the  wall  material.  Notice  that  Eqs.  (3.6)  and 
(3.10)  yield 

nf  =  (l-6c)ns.  (3-11) 

For  C  =  (r  —  rw)/Xci1  with  =  ^e0 Tp/e2riiQ1  the  natural  spatial  variable  in  the  sheath, 

U  (e<p)  =  npTp  +  ( ns  +  nf)mev^s  +  nirriiV 2ri  (3-12) 


the  Sagdeev’s  potential,  and  applying  the  boundary  condition 

d(e(j)/Tp ) 


d( 


Q 


0, 


(3.13) 
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the  integration  of  the  quasiplanar  Poisson  equation  yields 


(3.14) 

Two  other  boundary  conditions  at  point  Q  are  needed  to  close  the  sheath  problem.  These  are 
plasma  quasineutrality, 

0  =  U'Q  =  npQ  +nsQ  +  nfQ  -  niQ,  (3.15) 

and  the  Bohm  condition 


0  =  W 


"Tin 


ns  +  rif 
mev?s 


rij  \ 
rriiVri ' 


Solving  this  last  one  for  the  ion  velocity,  one  has 


(3.16) 


m  v2  _ _ ‘2e<t,WQTp'niQ _ 

nQ  2e<pWQnPQ  —  Tp(nSQ  +  uj-q) 


(3.17) 


It  will  be  shown  that  the  last  term  in  the  denominator  of  Eq.  (3.17)  is  small  always,  leading  to 
VriQ  -  y/Tp/rrii. 

One  magnitude  of  interest  is  the  heat  deposited  at  the  walls  by  the  currents  of  ions  and 
primary  electrons.  This  energy  flux  (per  unit  of  channel  length)  satisfies 

Qrw  =  2TpgrpQ  +  -mtv^.tWgnQ  =  ^  ^ — b  - mivliQ  +  e^wQ^griQ,  (3.18) 


(for  zero  axial  flow  of  the  plasma). 


3.3.1  Dimensionless  solution 

Parameters  Tp  and  riiQ  are  used  to  non-dimensionalize  the  sheath  equations;  hats  will  be  used 
to  call  non-dimensional  variables  (i.e.  <p  =  e(j>/Tp).  Then,  the  dimensionless  sheath  solution 
depends  on  just  two  parameters,  Sw  and  5C.  Figures  3.2(a)  and3.2(b)  show  the  sheath  potential 
profile  for  different  values  of  5C  and  Sw.  For  5W  given  [Fig.  3.2(a)],  the  net  SEE  flux, 


9rfQ  T  grsW  —  ^c^w9rpWj  (3.19) 

decreases  with  the  collision  parameter,  Sc.  Then,  in  order  to  satisfy  the  zero-current  condition, 
which  can  be  expressed  as 


9riQ  =  (1  -  5c5w)npQ  exp  6^WQ  yC^-,  (3.20) 

a  larger  (f>qW ,  is  required.  A  similar  reasoning  explains  the  evolution  of  the  profiles  with  Sw  in 
Fig.  3.2(b);  this  figure  can  be  compared  with  a  similar  one  for  Sc  =  1  in  Ch.  2. 

Solving  for  <j>wQ,  Eq.  (3.20)  yields 


^WQ 


e(t>WQ  _  jn 

T 


+  ln(l  -  5C5W )  +  In 


(3.21) 
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Then,  as  long  as  1  —  5C5W  <  0(y/'me/mi)1  one  has  4>wq  —  In  \/?7q/27rme.  A  very  interesting  fact 
is  that  4>wq  depends  on  5C  and  Sw  mainly  through  the  product  SCSW.  Physically,  it  means  that 
( t>wQ  is  fixed  by  the  ratio  grpQ/9riQ  =  1  —  SCSW.  Since,  in  addition,  the  typical  Debye  length  of 
the  sheath  does  not  depend  practically  on  the  SEE,  cases  with  the  same  value  of  SCSW  are  going 
to  have  very  similar  sheath  profiles.  Figures  3.3(a)  and  3.3(b),  show  density  profiles  for  two 
very  dissimilar  cases,  and  confirm  that  conclusion;  in  Fig.  3.3(b),  notice  that  the  SEE  density 
modifies  only  a  local  region  near  the  wall  and  has  no  global  impact  on  the  sheath  solution. 


3.3.2  Charge-saturation 


When  the  product  5C5W  becomes  close  to  one,  4>wq  decreases  and  the  sheath  approaches  the 
charge-saturation  limit  (CSL).  This  is  characterized  by  the  condition  d<p/d(\w  =  0  or,  using 
Eq.  (3.14),  by 

0  =  U(e(f>w)  ~  U{e(j)q).  (3.22) 

This  condition  yields  the  relation,  5W  =  <^(5C),  the  two  sheath  parameters  satisfy  at  the  CSL. 
Figure  3.4  shows  that  5*m  increases  as  5C  decreases,  and  attains  very  large  values  for  5C  <C  1;  one 
has  1  —  6^  ~  for  5C  =  1,  and  5 ^  ~  0 . 24 \Jmt  jrne  for  5C  =  0.  For  real  materials, 

the  maximum  SEE  yield,  5wm,  is  between  1  to  2  for  many  metals,  and  of  the  order  of  10  for 
some  dielectrics  [55,  43].  This  leads  to  the  important  conclusion  that  the  sheath  does  not  reach 
charge-saturated  conditions  in  low-collisional  plasmas,  more  precisely,  for  5C  <  l/5wm. 

For  wall  and  plasma  conditions  leading  to  5W  >  S*W(5C)  the  sheath  is  in  the  charge-saturated 
regime(CSR):  a  potential  well  is  formed  near  the  wall  to  reduce  the  actual  flux  ratio  of  popu¬ 
lation  s  reaching  the  bulk  of  the  plasma  to  5*w.  For  Tw/Tp  — *  0,  the  potential  well  is  negligible 
and  the  dimensionless  sheath  solution  for  the  CSR  is  the  same  that  for  the  CSL. 

The  CSL  relation  5W  =  S^(SC)  is  explained  from  the  pressure  balance  implicit  in  Eq.  (3.22). 
Extending  the  analysis  of  Ch.  2  to  Sc  <  1,  the  dominant  terms  in  Eq.  (3.22)  at  the  CSL  yield 


t rn  ^  \ga\  +9f  _  K (2  -  5C) 
V-e  9i  ~  1  -  <5C<5* 


(3.23) 


Indeed  the  approximation 

l  Z  yj,  =  2.35  +  5.644  (3,24) 

^ W\ ^  ^c) 

fits  very  well  the  exact  solution  for  5^(5C);  both  are  indistinguishable  in  Fig.  3.4.  Finally, 
Eq.  (3.23)  yields 

\gs I  ~  ^  >  9i  if  sc  =  0(1),  /3  25n 

|^|  »  gp  cz  gi  if  5C  <  1.  1  ' 

indicating  that,  at  the  CSL,  the  electrical  current  balance  among  the  different  species  depends 
on  the  SEE  collisionality. 


3.3.3  Results 


Figures  3.5(a)-3.5(c)  show  the  influence  of  5W  and  5C  on  some  important  sheath  magnitudes;  in 
Fig.3.5(c) 


/v  _  Qrw 

Qrw  = 


Tpg< 


riQ 


1  -  5W5C 


+  <t>WQ  + 


2  (1  +  2<PWq)  (2  Sc)nSQ 


(3.26) 
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is  the  (dimensionless)  energy  flux  deposited  at  the  wall  per  ion  impinging  it  and  unit  of  channel 
length.  The  region  of  valid  sheath  solutions  turns  out  to  be  bounded  by  four  limit  cases:  (i) 
total-trapping,  (ii)  zero-trapping,  (iii)  maximum  SEE  yield,  and  (iv)  sheath  charge-saturation. 

Figures  3.6(a)-3.6(c)  are  similar  to  the  preceding  ones  but  emphasize  that  the  role  of  the 
plasma  temperature  Tp  and  a  maximum  SEE  yield.  For  5W(TP )  we  used  the  Lye-Dekker  formula 
[referenced  by  Vaughan [55]] 

W)  =  ,  (3.27) 

J \zm) 

where  Tpm  is  the  electron  temperature  where  5wm  is  reached,  f(z )  =  2Ta35[l  —  exp(— z1'35)], 
zm  ~  1.84,  and  f(zm )  ~  0.725. 

3.4  PLANAR  PRESHEATH  ANALYSIS 

The  presheath  problem  is  formulated  and  solved  as  in  Ch.  2,  so  that  we  can  avoid  to  justify  every 
point.  The  plasma  in  the  presheath  consists  of  ions(i),  the  two  symmetrical  populations(/)  of 
free  SEE,  and  a  population  e  of  confined  electrons.  Th 


as  the  transition  (or  Bohm)  condition  at  the  presheath  side  of  point  Q.  Both  Eqs.  (3.17)  and 
(3.33)  are  generalized  forms  of  the  sonic  Bohm  condition  for  multi-species  plasmas.  Physically, 
the  sonic  condition  on  a  single  fluid  is  related  to  the  balance  between  convection  and  pressure. 
In  this  multi-species  plasma,  the  quasineutral  electrostatic  interaction  of  the  ions  with  several 
electron  populations  is  felt  like  a  pressure  term  by  the  ion  fluid,  the  right  sides  of  Eqs.  (3.17) 
and  (3.33)  being  the  effective  pressure  (divided  by  n^Q )  on  the  ion  fluid. 

Imposing  vriq  to  satisfy  the  right-hand  sides  of  Eqs.  (3.17)  and  (3.33),  the  relation  between 
Te  and  Tp  we  sought  is  obtained, 


neQ  _  npQ  _  nfQ  /g  g^\ 

Te  Tp  2e<pwQ 5 

The  conservation  of  the  electron  flux  and  density  at  point  Q, 

neQ  =  nPQ  +  5cnSQ,  geQ  =  gpq  +  ScgsQ,  (3.35) 

completes  the  matching  conditions  for  the  electron  populations  in  sheaths  and  presheath. 

Since  the  sheath  problem  is  solved  already,  it  is  convenient  to  integrate  the  presheath  prob¬ 
lem  from  point  Q  towards  point  M,  with 


7]  = 


(Tq  -  r)vu 

\/Tp/mi 


(3.36) 


as  independent  variable,  and  using  and  Tp  to  non-dimensionalize.  In  this  way,  apart  from 
5W  and  Sc.  the  presheath  solution  depends  on  one  extra  parameter,  ar  —  vr/vw. 

The  integration  ends  at  vriM  =  0  (or  d<j)/dr\M  =  0)  and  defines  the  value  of  g  at  point  M, 

Vm  Qy)  (3.37) 

Since  the  channel  width  hc  is  known,  gM  determines,  in  fact,  the  production  frequency, 

_  2gM\/Tp/mi 


hr 


(3.38) 


Using  the  characteristic  frequency  =  hc1  and  calling  Da  =  vOLjv{ 0,  ( a  —  w,r )  the 

functional  relation  (3.37)  becomes 


U/j  r :  ^vj ?  & c)  (3.39) 

which  is  the  plasma  balance  equationoi  this  plasma,  and  shows  that  the  net  production  frequency 
is  not  an  independent  parameter  [54],  For  total-thermalization  (5C  =  1)  vw  did  not  depend  on 
5W,  just  on  7y.  For  5C  <  1,  both  5W  and  5C  affect  vw.  However,  since  the  density  of  the  /- 
populations  is  small  [according  to  Fig.  3.5(b),  2 is  always  less  than  15%]  the  influence 
of  the  free  SEE  on  the  presheath  solution  is  small.  Figure  3.7  plots  vw{pr )  for  different  values 
of  (SW,SC),  showing  that  the  influence  of  these  two  parameters  is  rather  irrelevant. 
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3.5  CONCLUSIONS 

The  model  presented  here  assumes  that  if  SEE  is  partially-collisional,  it  can  be  split  in  two 
groups,  one  trapped  and  eventually  thermalized  with  the  main  plasma,  the  second  one  consti¬ 
tuted  by  two  free  counterstreaming  beams,  which,  for  the  planar  case,  are  collected  back  at  the 
wall  opposite  to  the  emission  one. 

The  density  of  the  two  free  SEE  beams  is  low  in  most  of  the  sheath  and  the  whole  presheath. 
Then,  the  presheath  solution  is  practically  identical  to  the  total-trapping  case.  The  ion  flux 
into  the  walls  is  determined  by  presheath  properties:  plasma  temperature,  channel  width,  and 
one  contribution  from  the  axial  behavior. 

Net  SEE  flux  increases  with  the  product  SC5W,  and  reduces  the  flux  of  primary  electrons, 
grpw-,  needed  to  keep  the  zero-current  condition.  The  sheath  potential  drop  is  self-adjusted 
to  collect  the  required  gTpw-  As  a  consequence,  the  sheath  solution  is  almost  invariant  for 
5W5C  =  const.  This  explains  that  charge  saturation  requires  a  larger  SEE  yield  as  the  SEE 
collisionality  in  the  plasma  decreases.  Indeed,  for  low  SEE  collisionality,  the  maximum  SEE 
yield  of  the  wall  material,  5wm,  will  be  below  5*m.  and  charge  saturation  cannot  be  reached. 

Finally,  the  collisionless  or  zero-thermalization  case  does  not  pose  any  particular  problem. 
For  all  practical  cases,  the  sheath  solution  is  the  same  than  for  the  zero  SEE  case. 
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(a) 


(b) 


presheath 


Figure  3.1:  Sketch  of  the  radial  models  in  (a)  presheath  and  (b)  sheath. 
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c 


c 


Figure  3.2:  Sheath  solution.  Prohles  of  the  electrostatic  potential  for:  (a)  5W  =  0.9  and  several 
values  of  Sc  from  0  tpilt 


(a)  1 


Figure  3.3:  Sheath  solution.  Profiles  of  the  different  density  species  for:  (a)  5W  =  0  and  Sc  =  1; 
(b)  5W  =  2  and  5C  =  0. 


5C 


Figure  3.4:  Sheath  solution.  The  charge-saturation  limit. 
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Figure  3.5:  Sheath  solution,  (a)  Potential  drop,  (b)  relative  free  SEE  density  of  at  the  sheath 
transition,  and  (c)  radial  energy  losses  in  terms  of  5W  and  Sc.  up  to  the  CSR. 
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Chapter  4 


AXIAL  MODEL  WITH  HEAT 
CONDUCTION 

4.1  INTRODUCTION 

Recently  we  presented  a  one-dimensional(lD)  model  (hereafter  called  Model  I)  of  the  stationary 
axial  structure  of  the  plasma  discharge  inside  a  Hall  thruster  and  in  the  near-plume[10].  The 
structure  consisted  in  an  anode  sheath,  a  region  of  free-diffusion  for  electrons  and  reverse  flow 
for  ions,  the  ionization  zone,  and  the  acceleration  region,  which  extends  into  the  plume,  Fig.  4.1. 
Although  quantitative  accuracy  could  not  be  expected  from  Model  I,  it  reproduced  rather  well 
most  of  the  main  features  of  the  experimental  data  [15,  31,  38,  32],  Compared  to  previous 
ID  models[39,  28,  19],  Model  I  unveiled  the  combined  importance  of  the  electron  pressure  and 
the  reverse  flow  of  ions  in  shaping  the  plasma  structure.  In  addition,  the  inclusion  of  the 
near-plume  (with  the  neutralization  surface)  allowed  closure  of  the  formulation  with  realistic 
boundary  conditions,  and  parametric  studies  on  thruster  performance.  Of  particular  interest 
was  the  determination  of  the  parametric  domain  of  stationary  solutions,  which  was  bounded 
by  the  case  of  zero  ion  reverse  flow  [6,  8]. 

The  main  defect  of  Model  I  came  from  the  simplified  treatment  of  the  electron  energy 
equation,  which  included  neither  heat  conduction  nor  the  losses  coming  from  the  interaction 
with  the  lateral  walls  of  the  chamber.  The  resulting  temperature  profile  had  too  large  gradients 
around  the  ionization  region  and  a  too  large  peak  temperature.  A  hybrid  (particle/fluid) 
numerical  model  by  Fife,  which  includes  heat  conduction  and  wall  losses,  yields  smoother 
temperature  profiles[25].  In  this  chapter  we  present  a  ID  macroscopic  model  (hereafter  called 
Model  II)  which  adds  heat  conduction  to  Model  I.  The  first  point  to  be  emphasized  is  that 
heat  conduction  is  not  a  new  forcing  term  added  to  the  equations  of  Model  I:  it  transforms  the 
electron  energy  equation  from  first  to  second  order  for  the  electron  temperature,  a  fact  having 
major  implications  in  the  mathematical  treatment  and  the  parametric  domain  of  stationary 
solutions.  The  effects  of  the  plasma  interaction  with  lateral  walls  are  not  included  in  Model 
II.  Two  are  the  reasons:  first,  a  detailed  radial  model  of  that  interaction,  to  be  coupled  to  the 
axial  model,  is  still  under  investigation [9];  and  second,  it  is  of  interest  to  understand  the  plasma 
response  in  thrusters  where  lateral  losses  are  small. 

The  chapter  includes  the  detailed  formulation  of  Model  II,  the  discussion  of  the  spatial 
response,  and  a  parametric  investigation  on  the  thruster  stationary  performances.  In  addition, 
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three  issues,  not  treated  in  detail  in  Ref. [10],  are  analyzed  here.  The  first  one,  treated  more  or 
less  explicitly  in  the  recent  literature,  is  the  evaluation  of  the  diffusive  approximation  for  the 
electron  motion.  Tahara  et  al.  [53]  and  Barral  et  al.  [13]  have  added  different  inertial  terms 
into  their  ID  macroscopic  models;  Fedotov  et  al.  [24]  suggest  the  presence  of  unmagnetized  and 
magnetized  populations  of  electrons  near  the  thruster  exit;  Haas  and  Gallimore[32]  conclude, 
from  their  experiments,  that  the  drift  energy  can  be  high  enough  to  contribute  to  ionization. 
The  second  issue  is  the  computation  of  the  ion  temperature,  seeking  a  confirmation  that  it  can 
be  neglected  everywhere  compared  to  either  the  electron  temperature  or  the  ion  beam  velocity. 
And  the  third  one  is  a  detailed  evaluation  of  the  thrust:  computations  based  in  the  momentum 
thrust  of  the  plasma  either  at  the  channel  exhaust  (as  we  did  in  Ref.  [10]),  or  in  the  far-plume 
(as  simulation  models  usually  do[25])  are  not  accurate  enough. 

The  chapter  is  organized  as  follows.  In  Sec.  II  we  formulate  the  quasineutral  model, 
including  a  detailed  discussion  on  the  diffusive  approximation  for  electrons.  In  Sec.  Ill  we 
discuss  the  model  singular  points,  the  anode  sheath,  the  boundary  conditions,  and  the  numerical 
procedure.  In  Sec.  IV  we  analyze  the  spatial  solution,  the  different  collision  frequencies,  the 
validity  of  the  diffusive  approximation,  and  the  ion  temperature.  In  Sec.  V  we  discuss  the 
thruster  performance  in  the  parametric  domain  of  stationary  solutions.  Conclusions  are  written 
in  Sec.  VI. 


4.2  THE  QUASINEUTRAL  MODEL 


Geometrical  sketches  of  the  thruster  and  the  axial  model  are  drawn  in  Fig.  4.1(a)-4.1(b). 
The  general  hypotheses  of  the  present  Model  II  are  the  same  as  in  Model  I.  In  the  ID  axial 
approximation,  plasma  variables  represent  average  values  on  each  cross-section  and  depend 
only  on  the  axial  variable  x.  The  macroscopic  formulation  will  consist  of  fluid-like  equations 
for  ions(i),  electrons(e)  and  neutrals(n).  The  thruster  channel  is  of  length  L  and  radial  area 
Ac.  The  plasma  discharge  extends  outside  the  thruster  into  a  plasma  plume  of  radial  area 
A(x)  to  be  determined.  The  magnetic  field  is  considered  purely  radial,  B  =  B(x) lr,  with  the 
maximum  field  placed  near  the  thruster  exit,  and  two  different  semi-Gaussian  axial  profiles  will 
be  used  inside  the  thruster  and  in  the  near-plume.  Electrons  are  injected  into  the  plume  at 
a  neutralization  surface  (point  P)  placed  at  a  distance  LEp  from  the  channel  exit  (point  E). 
The  voltage  difference  between  anode  (point  A)  and  point  P  is  the  discharge  voltage,  V^,  and 
the  electron  current  delivered  at  the  neutralization  surface  is  the  discharge  current,  Id-  One 
part  of  this  current  flows  outwards  and  neutralizes  the  ion  current  1^ ;  the  other  part,  Id  —  Iip1 
diffuses  inwards  across  the  radial  lines  of  the  applied  magnetic  field  and  ionizes  the  mass  flow  of 
neutrals,  m^,  injected  at  the  anode;  subscript  oo  will  refer  to  downstream  conditions  far  away 
from  the  cathode  P.  No  effects  of  the  radial  interaction  of  the  plasma  with  the  lateral  walls  of 
the  chamber  will  be  included  in  this  model. 

The  axial  plasma  flow  is  considered  quasineutral  everywhere  except  in  a  thin,  electron- 
repelling  sheath  attached  to  the  anode  (region  AB  in  Fig.  4.1,  with  xp  —  %a  =  0  in  the 
quasineutral  scale).  The  potential  jump  in  the  sheath,  4>ab  =  4>b  ~  4>a  >  0,  adjusts  the  small 
diffusive  electron  flow  in  the  channel  to  the  thermal  flow  collected  at  the  anode.  Since  all  ions 
are  created  by  ionization  within  the  thruster,  plasma  quasineutrality  implies  back-flow  of  ions 
in  the  rear  part  of  the  channel. 
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4.2.1  Electron  dynamics 

We  review  here  the  closed-drift,  diffusive-motion  model  in  order  to  evaluate  the  importance 
of  inertial  effects  in  the  electron  dynamics.  Also,  the  diffusive  equation  for  heat  conduction  is 
formulated. 

The  vector  equation  for  the  electron  momentum  is 

meneve  ■  Vue  =  —  VneTe  +  eneV<fi  —  eneve  A  B  —  Re.  (4.1) 

Apart  from  conventional  symbols, 

Re  =  mene(veve  -  uelvt  -  venvn ) 
is  the  friction  force  due  to  collisions,  with 

—  ven  +  lSei  +  van01  (4-  2) 

the  total  collision  frequency  for  electrons.  This  includes  contributions  from  electron-neutral  col¬ 
lisions  (ven),  electron-ion  collisions  (vei),  and  anomalous  effective  collisions  representing  Bohm- 
like  diffusion  [20],  which  is  expressed  as 


^ ano 


with  u oe  =  eB/me  the  electron  gyrofrequency  and  a;#  an  empirical  parameter  (with  a  ’classical’ 
value  of  1/16).  When  current  losses  to  lateral  walls  are  taken  into  account,  ve  has  another 
contribution  from  near- wall-conductivity  [21,  9]. 

In  the  1-D  axial  approximation,  the  electron  velocity  held  is  ve  =  vxelx  +  vge lg,  and  spatial 
gradients  reduce  to  the  axial  direction,  V  =  1  Xd/dx.  Then,  applying  the  closed-drift  hypothesis, 


ue  >  ue,  (4.3) 

and  the  ansatz  \vge/vxe\  ~  ue/^e,  Eq.  (4.1)  simplifies  into 

-^-neTe  +  ene^,  (4.4) 

Cvtv  Cvw 

UeVxe  -  VeVge.  (4.5) 

The  hrst  term  in  Eq.  (4.5)  is  the  inertia  of  the  electron  azimuthal  motion  and  is  the  dominant 
convective  contribution.  Then,  the  diffusive  approximation  requires 


0  ~  meneuievge  - 
dvee 

vxe—r~  -  - 
ax 


1 

dVge 

U)e 

dx 

<  1, 


and  reduces  Eq.  (4.5)  to 

Vde/Vxe  -  -Ue/Ve. 

The  equation  for  the  electron  internal  energy  reads 
'3 


V  •  (-■ neTeve  +  =  -neTeV  •  ve  -  ve  ■  Re  -  h>ineaiEi 


(4.6) 

(4.7) 


(4.8) 
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Here,  qe  is  the  electron  heat  conduction, 


Ve  ■  Re  —  menevev9e 

is  the  Joule  heating  (we  assumed  vn  <C  vXi  ~  vxe  <C  vge ) ,  and  the  last  term  on  the  right  side 
accounts  for  ionization  and  excitation  losses;  £)  is  the  ionization  energy,  cq  is  the  ionization 
cost  factor  (taking  into  account  excitation  and  multistep  ionization),  and  i \  =  nnRi(Te )  is  the 
ionization  frequency.  Only  single  ionization  will  be  considered,  Dugan’s  formula[22] 

„  ,  12  Ei 

<*i(Te)  -  2+  4  exp— , 

will  be  used,  and  the  expression  for  the  ionization  rate,  Ri(Te),  will  be  taken  from  Ref.[10]. 

Equation  (4.8)  for  the  internal  energy  identifies  clearly  the  effects  contributing  to  the  electron 
heating.  Thus,  it  is  more  adequate  to  use  than  the  equation  for  the  electron  total  energy,  where 
the  contributions  to  the  internal  and  mechanical  energy  are  mixed.  Indeed,  the  equation  for 
the  mechanical  energy  comes  out  from  Eqs.  (4.4)-(4.5): 

d  (l  2  \  d  rp  2  # 

neVxe—  ^-meveeJ  =  -vxe—neTe  -  vemenev9e  +  enevxe—.  (4.9) 


-neTe  -  vemenev9e  +  enevx 


Notice  that  the  term  on  the  left  side  corresponds  to  the  convective  contribution,  and  must  be 
neglected  when  the  diffusive  limit  (4.7)  is  used. 

Tahara  et  al.  [53]  incur  an  inconsistency  by  including  the  electron  mechanical  energy  in 
the  total  energy  equation  while  using,  at  the  same  time,  the  diffusive  approximation  in  the 
momentum  equation.  Barral  et  al.  [13]  proposed  a  non-stationary  model  which,  in  the  stationary 
limit,  reduces  to  Eqs.  (4.4)-(4.5)  and  (4.8).  However,  (i)  in  the  non-stationary  formulation, 
they  omit  the  temporal  term  dv9e/dt  associated  to  the  azimuthal  inertia,  and  (ii),  it  is  unclear 
whether  they  really  include  the  inertia  term  in  the  numerical  integration,  since  they  omit 
explicitly  a  boundary  condition  for  v9e. 

To  close  the  equations  for  electron  dynamics,  a  transport  equation  is  needed  for  qe.  In  the 
diffusive  approximation,  that  equation  is[17] 


-neTe  VTe  +  eqe  A  B  +  meveqe  a 

from  where  the  azimuthal  and  axial  components  of  qe  satisfy 

0  ~  -ueqxe  -  ueqde, 

n  5 Pe  dTe  , 


2  mP  dx 


+  ueqge. 


(4.10) 

(4.11) 


The  resulting  equation  for  qxe  coincides  with  the  one  obtained  by  Fife  [25]  from  the  postulate 
of  equal  diffusitivities  for  heat  and  mass  transport. 


4.2.2  Quasineutral  equations 

Equations  for  ions  and  neutrals  are  the  same  as  in  Model  I.  Adding  them  to  the  above  diffusive 
model  for  the  electrons,  the  stationary  model  for  the  quasineutral  plasma  consists  of 


-(Anevxi)  =  ^-^(Anevxe)  =  -  j^(An„u„) 


(4,12) 
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vn  =  const, 

Id..  9  ,  d(j) 

—  —  ( Aminevxi )  =  -ene  —  +  minei>ivr 

n  d  rp  .  d(P 

0  =  — —neTe  +  ene- - vdrnenevxe. 

dx  dx 


dneTj 
dx  ’ 


(l 

Adx  \2 
dTe  _  2 mevd 

dx  5neTe 


dvT 


(  7;TenevXe  +  qxe  )  =  -neTe^AA  +  vdmenev2xe  -  VinectiEh 


Qxe 


(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 


Symbols  are  conventional  or  were  defined  already  in  Ref.  [10].  Equations  (4.7)  and  (4.10) 
have  been  used  to  eliminate  vge  and  qge  in  the  Ohm  law  (4.15)  and  the  Fourier  law  (4.17), 
respectively.  The  axial  discharge  is  governed  by  two  effective  frequencies:  the  axial  diffusion 
frequency  for  the  magnetized  electrons, 


^  =  we2M,  (4.18) 

and  the  ionization  frequency  =  nnRi(Te).  The  area  A(x)  and  thickness  h(x)  of  the  annular 
cross-section  of  the  plasma  jet,  remain  constant  within  the  channel,  and  they  satisfy 

—  In  .4  =  —  lnh  =  —  tan  5  =  — (4.19) 

dx  dx  h  vxih 

in  the  plume;  here  5(x)  is  the  ’mean’  semi-angle  of  divergence  of  the  plume  and  cpi  is  the  velocity 
of  lateral  expansion  of  the  plume,  which  we  will  define  below. 

To  close  the  set  of  Eqs.  (4. 12)- (4.19),  an  equation  for  the  ion  temperature  is  needed.  For  the 
moment,  we  apply  the  ansatz  Tt  <g;  Te,  which  allows  us  to  neglect  the  ion  pressure  in  Eq.  (4.14) 
and  thus  to  drop  Tt  from  the  equations.  The  ion  temperature  and  the  validity  of  the  ansatz 
will  be  discussed  later. 

The  quasineutral  model  applies  between  the  entrance  to  the  anode  sheath  (point  B)  and 
x  =  +oo,  except  for  a  discontinuity  on  the  electron  current  at  the  neutralization  surface  (point 
P).  The  anode  sheath  must  be  solved  separately,  in  its  own  natural  scale  (the  Debye  length), 
and  provides  boundary  conditions  for  the  quasineutral  model  at  point  B. 

Some  magnitudes  of  interest  are  the  particle  flows  TQ  =  Anava,  (a  =  i,  e,  n);  the  anode 
mass  flow,  ttia  =  rnjTm,  and  the  discharge  current,  Id  =  eTd,  with 

T  m  =  Ti(x)  +  Tn(x)  =  const, 

rd  =  Tj(x)  -  Te(x)  =  const; 
the  utilization  efficiency  and  the  ionization  fraction, 

Vu  =  ?7ioo,  r}i(x)  =  Ti(x)/Tm\ 


and  the  current  parameter, 


mjd  r  d 


4.3 


MODEL  INTEGRATION 


4.3.1  Singular/sonic  points 


Solving  Eqs.  (4.12)-(4.19),  with  T)  =  0,  for  the  spatial  derivatives  we  obtain  a  matrix  relation 
of  the  form 


dY 

a  -m2)— =  F(n 


dx 


(4.20) 


where 


M 


y/Tjm 


is  the  isothermal  Mach  number  for  the  ion  axial  flow,  Y  groups  the  eight  plasma  variables 
ne,  nn,  vXi,  vxe,  vn,  Te,  qxe.  and  <p  ,  and  F(Y)  is  a  regular  function.  In  particular,  the  scalar 
equation  for  vXi  can  be  written  as 


dV7 


dx 


Vi  —  vr 
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-vdmevxe(l 


5  neTev1 
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Te(  1  -  M2)  +A~dx. 


2vxi  -  v„)  +  rriiv: 


din  A 


2  _ 

xi  dx 


(4.21) 


Sonic  points,  M  =  ±1,  are  singular  points  of  Eq.  (4.20).  Notice  that  the  definition  of  the  Mach 
number  in  Model  II  does  not  include  the  specific  heat  ratio  5/3  of  Model  I,  and  function  G  is 
different  too.  As  in  Ref.  [10]  we  will  look  for  solutions  with  a  regular  sonic  transition, 

Gs  =  0  and  Ms  =  1  (4.22) 

inside  the  thruster  (point  S  in  Fig.  4.1)  and  a  supersonic  plasma  flow  at  the  thruster  exhaust. 


4.3.2  Anode  sheath 

The  Debye  sheath  attached  to  the  anode  (region  AB  in  Fig.  4.1)  completes  the  1-D  model  of 
the  channel  and  defines  boundary  conditions  at  point  B  for  the  quasineutral  equations.  Under 
usual  conditions,  the  sheath  is  electron-repelling  and  the  space-charge  held  adjusts  the  potential 
jump  in  the  sheath,  <f> ab  =  4>b  ~  4>a  >  0,  to  a  value  such  that  the  how  of  electrons  reaching  the 
anode  is  equal  to  the  diffusive  how  coming  from  the  quasineutral  channel. 

In  the  distinguished  limit  A D/L  — >  0,  the  problem  in  the  sheath  reduces  to  collisionless, 
conservation  equations.  In  particular,  particle  and  energy  hows  are  constant.  Assuming  a 
quasi-Maxwellian  electron  distribution  with  temperature  TeB,  electron  magnitudes  at  points  A 
and  B  are  related  by 


(jleVXe)  B  — 


1 

--neiJ  exp 


—e(f>AB 

TeB 


{  2nevxeTe  T  Qxe) B  —  (2 TeB  +  &(j) Ab)  (Ae'VXe) B • 


(4.23) 
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The  matching  of  this  ion-attracting  sheath  with  the  quasineutral  plasma  at  point  B  requires 
the  plasma  to  verify  the  Bohm  condition 


Mb  =  -1,  (4.24) 

which  states  that  the  ions  must  flow  into  the  sheath  with  the  plasma  local  sound  speed.  Using 
Eqs.  (4.23)  the  potential  jump  across  the  sheath  satisfies 


e4>AB 


=  In 


-  eB 


'  eB 


2^meV2xeB 


=  1„( 


rrij 
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2tt  me  1  T. 


eB 


(4.25) 


and  the  electron  energy  flux  deposited  by  conduction  into  the  anode  sheath  is 


QxeB  —  (jl'e'VxeTe)  B  ( 


e(f>AB  1 


'  eB 


)■ 


(4,26) 


4.3.3  Boundary  conditions 

Following  Refs.  [10]  and  [44]  the  appropriate  velocity  for  the  plume  radial  expansion  is  the 
plasma  sound  velocity  at  the  thruster  exhaust, 


cpi  =  y/TeE/mi.  (4.27) 

Equations  (4.12)-(4.17)  need  eight  boundary  conditions.  Seven  of  them  are  similar  to  Model 
I  and  there  is  one  additional  condition  for  qxe.  The  eight  conditions  are 

i)-ii)  The  injected  flow  and  the  velocity  of  neutrals  at  the  anode,  m,A  and  vnA(=  vnB),  are 
known. 

iii)  The  electron  temperature  at  the  neutralization  surface,  TeP,  is  known. 

iv)  The  discharge  voltage,  =  4>a~  (t)p-  is  known. 

v)  There  is  a  regular  sonic  point  inside  the  channel  characterized  by  Eq.(4.22). 

vi) -viii)  The  matching  with  the  anode  sheath  yields  conditions  (4.24)-(4.26). 

4.3.4  Integration  procedure 

From  the  preceding  equations  and  conditions,  it  turns  out  that  the  quasineutral  model  requires 
input  parameters  of  different  kinds.  There  are  design  parameters  (Ac.  L).  control  parameters 
(Vd,  B{x ),  ttia ),  plasma  parameters(Tep,  vnB ),  and  ’empirical’  parameters  ( LEp ,  aB).  The 
main  output  parameters  are  the  discharge  current  /<*,  the  thrust  F ,  and  the  thrust  efhciency 
77,  the  two  last  ones  to  be  defined  later.  Other  output  parameters  of  interest  are  the  position 
of  the  sonic  point,  xB,  the  divergence  angle  at  the  thruster  exhaust,  SE,  and  several  partial 
efficiencies. 

Plasma  equations  are  non-dimensionalized  as  in  Model  I;  hereafter,  an  upper  tilde  will 
indicate  dimensionless  variables.  Then,  they  are  integrated  with  Runge-Kutta  plus  shooting- 
routines  from  point  S  towards  points  B  and  P  independently.  This  means  that  plasma  param¬ 
eters  at  point  S  are  used  as  input  parameters  of  integration.  A  Taylor  expansion  solves  the 
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indeterminacy  Gs/(  1  —  AfJ)  =  0/0  and  yields  the  slopes  of  plasma  variables  at  point  S.  Using 
Eq.  (4.21),  the  regularity  condition  (4.22)  can  be  expressed  as 

t  ^xe  2(y  xe  /  (5neTe) 

udme  2vxi  -  vn 

It  is  easily  seen  that  this  equation  relates  r]iS  to  TeS,  qxes,  and  id.  The  integration  uses  these  three 
quantities  as  input  parameters  instead  of  the  natural  ones,  TeP,  qxeB ,  and  Vd.  As  a  consequence, 
an  iteration  is  needed  to  reproduce  the  desired  operating  conditions.  An  additional  difficulty 
of  Model  II,  compared  to  Model  I,  is  that  the  solutions  departing  from  point  S  do  not  converge 
easily  to  the  boundary  conditions  at  point  B.  Thus,  we  were  forced  to  launch  solutions  from 
point  S  and  B  separately  and  to  match  them  at  an  intermediate  point. 

Once  the  region  between  points  B  and  P  is  solved  the  solution  is  continued  to  x  =  +oo  by 
just  exchanging  the  condition  TxeP-  =  TxiP  -  Fd  by  TxeP+  =  TxiP. 

4.4  ANALYSIS  OF  THE  SPATIAL  SOLUTION 

4.4.1  Axial  structure  of  the  discharge 

Figure  4.2  shows  the  axial  profiles  of  main  plasma  variables  for  an  SPT-100  type  of  thruster; 
design  and  control  parameters  are  given  in  the  figure  caption.  The  discharge  presents  a  similar 
structure  to  the  one  obtained  in  Model  I  and  sketched  in  Fig.  4.1,  consisting  of  the  anode 
sheath,  the  diffusion  region,  the  ionization  layer,  the  internal  acceleration  region,  and  the 
plume.  However,  heat  conduction  introduces  some  differences  of  relevance.  First,  there  is  the 
desired  smoothing  of  the  temperature  profile,  with  a  peak  temperature  around  60  eV,  instead  of 
the  90  eV  of  Model  I.  Second,  the  sharp  transition  between  the  diffusion  and  ionization  regions 
found  in  Model  I  has  been  smoothed  too,  yielding  a  shorter  diffusion  region  here.  Notice  that  a 
simple  dimensional  analysis  of  Eqs.  (4.15)-(4.17)  suggests  that  heat  conduction  is  of  the  same 
order  of  the  convective  energy  flow.  Figure  4.2(f)  shows  that  heat  conduction  is  even  larger 
than  the  convective  flow  in  certain  regions  (besides,  heat  conduction  is  towards  the  anode  in 
most  of  the  channel,  but  away  from  it  in  the  plume). 

Plasma  equations  and  boundary  conditions  seem  to  force  an  almost  monotonic  profile  of  qxe 
between  B  and  S.  Therefore,  qxe  at  point  S  is  negative  and  the  point  of  maximum  temperature 
(i.e.,  qxe  =  0)  is  placed  downstream  of  point  S,  near  the  channel  exit.  Function  G  of  Model 
II,  Eq.  (4.21),  shows  that  the  energy  balance  allowing  the  regular  sonic  transition  at  point  S 
involves  losses  by  ionization,  Joule  heating,  and  heat  conduction.  However,  it  turns  out  that 
point  S  is  practically  in  the  acceleration  region  and  the  ionization  contribution  to  the  above 
balance  is  small,  which  makes  a  difference  with  Model  I.  For  the  case  of  Fig.  4.2,  the  ionization 
layer  ends  at  xB  ~  15.3  mm  approximately,  and  we  have  x$  —  16.4  mm,  and  qxe  =  0  at 
x  ~  22.6  mm. 

Figure  4.3(a)  compares  the  contributions  of  different  collision  terms  to  the  electron  collision 
frequency  ve.  The  main  conclusion  is  that  Bohm-like  diffusion  provides  the  main  contribution 
downstream  of  the  ionization  layer,  a  result  supported  experimentally  by  ve  by  Meezan  and 
Capelli[42],  However,  we  still  adjust  aB  ’empirically’  to  obtain  a  satisfactory  solution;  the  value 
we  use  is  of  the  order  of  the  one  used  by  Fife[25],  and  both  are  much  smaller  than  the  ’classical’ 


(4.28) 
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value  of  Bohm.  Figure  4.3(a)  shows  that  e  —  i  collisions  dominates  over  e  —  n  collisions  in  the 
rear  part  of  the  channel;  this  is  due  to  the  high  electron  density  and  low  electron  temperature 
arising  when  the  ionization  region  is  thin. 

Figure  4.3(b)  depicts  the  two  frequencies  governing  the  axial  discharge,  and  ud  =  udme/rrii 
[with  the  mass  ratio  in  ud  coming  out  from  the  dimensional  analysis  of  Eqs.  (4.12)-(4.17)].  The 
plots  show  that  z^  <C  vd  even  at  the  peak  of  z'*.  This  makes  the  structure  of  the  diffusion  plus 
ionization  region  (region  BH  in  Fig.  4.1)  not  simple  to  analyze:  on  the  one  hand,  rji,  vxi ,  and 
the  total  length  Lah  of  the  region  are  governed  by  z^;  on  the  other  hand,  Eqs.  (4.15)  and  (4.17) 
yield  that  the  profile  of  ne  and  Te  are  controlled  mainly  by  vd.  Changes  of  (f>  in  region  BH  are 
small,  of  the  order  of  Te  there;  for  the  case  of  Fig.  4.2,  one  has  that  <f>  increase  by  1  V  from 
point  B  to  the  point  with  vxi  =  0,  and  the  voltage  drop  between  point  B  and  point  S  (already 
in  the  acceleration  region)  is  4>bs  —  23.1V  only. 

The  structure  of  the  internal  and  external  acceleration  regions  (region  HP  in  Fig.  4.1) 
is  simpler  because  ionization  is  negligible  there.  The  quasineutral  model  simplifies  to  the 
conservation  relations 
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(4.30) 

(4.31) 


plus  Eqs.  (4.15)  and  (4.17)  for  electron  momentum  and  heat  conduction.  These  two  last 
equations  show  that  f  (— vxe)meuddx  is  the  relevant  integration  variable  for  the  acceleration 
regions.  Quasineutrality  and  Eq.  (4.29)  state  the  electrons  must  counterflow  against  the  ions 
with  a  velocity 

~vxe  =  ( id/Vu  ~  1  )vxi.  (4.32) 

Adding  to  these  facts  the  conditions  that  the  ion  low  is  mostly  supersonic  and  most  of  the 
voltage  drop  takes  place  in  the  region  HP,  the  integration  of  Eq.  (4.15)  yields 

2 eVd  fP  (  id  A  meul 

m  JH  'Vu  '  rriiVe 

This  expression  suggests 

Vd/2  oc  (id  -  1  )^Lhp 

V e 

as  an  important  scaling  law  among  the  main  thruster  parameters  for  efficient  operation. 

Kim  derived  a  similar  law[38].  Then,  assuming  that  (i)  LHP  ~  const  (for  a  given  geometry), 

1  /2 

(ii)  id  —  1  ~  const,  and  (iii)  ve  ~  Vd  jhc  (from  an  estimate  based  in  the  dominance  of  wall 
collisions  and  hc  being  the  channel  width),  he  concluded  that  optimum  operation  required 
Bmax  oc  V^1'2 .  Curiously,  when  Bohm  diffusion  dominates  one  has  ue  ~  aBBmoa;  and  the  scaling 
law  Bmax  oc  Vd^'2  still  holds,  provided  that  the  condition  id  —  1  ~  const  applies.  We  will  come 
back  to  this  discussion  on  optimum  operation  in  Sec. 4. 6. 2. 


-dx. 
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4.4.2  Inertial  effects  on  electrons 


We  evaluate  now  the  assumptions  of  closed-drift  and  diffusive  electron  dynamics.  Figure  4.4 
plots,  for  the  case  of  Fig.4.2,  three  parameters  related  to  electron  dynamics:  the  inverse  Hall 
parameter,  the  drift-to-internal  energy  ratio,  and  the  convective-to-magnetic  time  ratio.  We 
see,  first,  that  the  closed-drift  approximation,  Eq.  (4.3),  is  verified  in  the  whole  domain,  even 
near  the  cathode,  where  the  magnetic  held  is  vanishing.  However,  the  diffusive  approximation, 
Eq.  (4.6),  fails  in  the  vicinities  of  anode  and  cathode.  Electron  convective  effects  around  the 
anode  come  from  vxe  oc  vxi  and  the  acceleration  of  the  plasma  having  to  meet  the  Bohm 
condition  at  point  B.  These  effects  are  rather  local  and  we  believe  they  are  not  significant 
globally.  On  the  contrary,  the  inertia  effects  in  the  plume  could  have  more  impact  in  the 
plasma  discharge  since  they  are  less  localized  and  the  magnetic  held  becomes  residual  there. 
From  Eq.  (4.5),  one  has 

), 

with  dvge/dx  <  0  in  the  plume,  yielding  thus  a  larger  value  of  vxe/vge  than  the  diffusive 
approximation.  This  can  be  interpreted  as  an  extra  collisionality  and  could  be  an  important 
and  independent  contribution  to  the  overall  ’anomalous’  diffusion.  However,  without  solving 
the  coupled  equations  (4.4)-(4.5)  together  with  the  other  ones,  it  is  difficult  to  predict  the 
overall  effect  of  the  convective  term  on  the  discharge  characteristics. 

In  spite  of  the  large  values  reached  by  the  convective  azimuthal  term,  Fig.  4.4  shows  that 
the  drift  energy  remains  small  in  the  whole  domain.  The  explanation  lies  in  the  azimuthal 
inertia  term  of  Eq.  (4.5)  being  balanced  by  the  azimuthal  component  of  the  Lorentz  force, 
which  yields  no  contribution  to  the  drift  energy.  In  any  case,  we  must  be  aware  that  vge  is  very 
sensitive  to  the  local  values  of  the  magnetic  and  electric  fields,  and  relatively  small  variations 
of  these  magnitudes  could  make  me  »2,/2  T.  1  locally,  as  it  has  been  found  in  the  experiments 
of  Haas  and  Gallimore[32],  Were  this  the  case,  the  drift  energy  should  be  taken  into  account 
to  compute  the  ionization  rate,  Ri. 

4.4.3  Ion  temperature 

The  evolution  equation  for  the  ion  temperature  is 

Vxi^x^Ti  =  ~~itxAVxi + Vi  [\mi(Vxi  ~ Vxn )2  +  \^Tn  ~  ’  (32)  (4-36) 

where  Tn  is  the  temperature  of  neutrals  (which  we  are  ignoring).  Equation  (32)  requires  a 
boundary  condition  for  Tt  and,  at  the  same  time,  introduces  a  new  singular  point  at  vxi  =  0. 
The  correct  boundary  condition  for  a  regular  solution  comes  from  making  zero  the  right-hand 
side  of  Eq.  (32)  at  vxi  =  0:  this  yields  Tt  =  (3Tn  +  rriiVxn) /5.  Again  a  Taylor-expansion  problem 
must  be  solved  to  determine  the  value  of  dTjdx  for  a  regular  crossing. 

Figure  4.5  plots  Tt(x).  obtained  from  the  integration  of  Eq.  (4.36)  using  vxl  {x)  and  A(x)  from 
the  cold-ion  solution  of  Fig.  4.2,  and  Tn/Te  — >  0.  It  confirms  that  T)  is  totally  negligible  and  the 
ion  beam  behaves  as  monoenergetic.  Clearly,  the  small  value  of  T)  and  the  drawback  of  a  new 
singular  point  at  vxl  =  0  makes  inadvisable  to  add  Eq.  (4.36)  to  the  set  of  Eqs.  (4.12)-(4.19). 
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The  low  ion  temperature  is  due  to  the  short  ionization  region  and  the  absence  of  lateral 
recombination.  The  experimental  Retarding  Potential  Analyzer  data  show  that  T)  ~  2  —  3  eV 
at  exit.  The  discrepancy  might  be  related  to  the  spreading  of  the  ionization  zone  when  there 
is  wall  recombination;  ion-ion  and  ion-neutral  collisions  could  contribute  too. 


4.5  THRUST  AND  EFFICIENCY 


4.5.1  Thrust 


The  electrostatic  interaction  between  ions  and  electrons  plus  the  magnetic  forces  on  the  electrons 
are  the  mechanisms  by  which  plasma  thrust  is  transferred  to  the  solids  parts  of  the  engine.  Since 
the  magnetic  held  extends  some  distance  outside  the  chamber  and  the  plasma  plume  presents 
some  radial  divergence,  neither  the  momentum  thrust  of  the  plasma  at  the  chamber  exhaust 
(point  E)  nor  in  the  far  plume  (point  oo )  evaluate  accurately  the  thrust. 

Assuming  the  engine  to  be  at  rest  by  some  thrust  balance  device  and  zero  external  pressure, 
the  axial  momentum  balance  for  the  system  constituted  by  the  thruster  plus  the  internal  plasma 
yields 

F  =  FpE  +  Fext ,  (4.37) 

where  F  is  the  thrust  (exerted  on  the  thruster  balance), 

Fp(x)  =  ( rriiriiV. 2  +  rriinnvl  +  neT^j  A  (4.38) 

is  the  total  axial  momentum  of  the  plasma  at  a  given  axial  section,  and 

/•OO  /*OC 

Fext  =  /  {—eneve)  A  B  ■  1  xA(x)dx  ~  j  enevgeB  A(x)dx  (4.39) 

J  E  J  E 

is  the  external  thrust  due  to  the  magnetic  force  on  the  electrons  in  the  plume.  Therefore,  the 
leakage  of  the  magnetic  held  outside  the  thruster  channel  yields  a  contribution  to  the  thrust, 
which  is  a  unique  characteristic  of  this  type  of  electric  thrusters. 

From  Eq.  (4.14)-(4.15),  the  axial  momentum  of  the  whole  plasma  satisfies 


dFp 

dx 


rpdA  ,  UA 

neTe—  +  eneVoeBA, 
dx 


(4.40) 


where  we  wrote  evgeB  ~  —vxeVd  for  the  magnetic  force  exerted  on  an  electron.  Integrating  Eq. 
(4.40)  along  the  plume  and  using  Eq.  (4.37)  one  has 


where 


F  =  F  -  D 

±  x  poo 


D, 


ext 


no 

JE 


ext i 

(4.41) 

dA 

(4.42) 

——dx 

dx 

is  a  drag  contribution  due  to  the  plume  lateral  expansion.  Expressions  (4.37)  and  (4.41)  yield 
the  difference  between  the  actual  thrust  and  the  momentum  thrust  at  points  E  and  oo.  One 
has 


F  —  F  w 

*  poo  ±  ph 


Fext  +  D, 


exti 


(4.43) 
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which  means  that  the  external  acceleration  of  the  ions  is  transferred  only  partially  to  thrust  (via 
ion-electron  electrostatic  interaction  plus  electron-thruster  magnetic  interaction).  A  relevant 
point  is  that  were  the  magnetic  field  zero  outside  the  thruster,  all  external  acceleration  would 
go  to  plume  radial  expansion.  For  the  case  of  Fig.  4.2  we  have  that  the  contribution  of  the 
external  magnetic  held  amounts  to  Fext  j F  ~  15.8%  whereas  the  external  acceleration  of  ions 
is  ( Fpoo  -  FpE)/F  ~  19.7%. 


4.5.2  Thrust  efficiency 

The  thrust  efficiency  is  defined  as 


V  = 


2 mAIdVd  ’ 


and  is  conveniently  factorized  into  four  terms 

V  =  VuVdVeVo, 

where  rju  is  the  utilization  efficiency,  and 

Iioc  _  Vu  miVxioo 

Vd  =  -j-  =  —  Ve=  Q  T/  , 

h  i>d  2eVd 


Vo  = 


(rriirivxi) 


2  5 
oo 


(4.44) 

(4.45) 

(4.46) 


are  efficiencies  for  discharge,  beam-energy,  and  other  effects  (plume  expansion,  anode  sheath,...), 
respectively.  For  the  case  of  Fig.  4.2,  one  has  rju  —  100%,  7]d  —  85.6%,  r/e  —  101.4%,  r/0  ~  93.4%, 
and  T)  —  81.1%;  the  thin  ionization  region  and  the  fact  that  the  maximum  voltage  is  not  at 
point  A  but  at  vxi  =  0  are  responsible  of  r]e  >  100%.  Notice  that  the  definition  of  r]  is  based  in 
the  anode  mass  how  and  does  not  include  the  power  spent  in  the  magnetic  circuit. 

The  discharge  efficiency  is  the  main  factor  reducing  the  thrust  efficiency.  The  relation 
between  (or  id )  and  the  ionization  losses  is  deduced  from  the  energy  balance  for  the  whole 
plasma.  From  the  model  equations,  the  energy  balance  for  the  plasma  is 

(4.47) 


IdVd  ~  ^  (Tqoo  —  WaA)  +  Qi 


a=i,e 


with 


1  5 

f  1  n  [^^(^^aba  ^Tv)  T  Q'sa] -4 


(4.48) 


VjOijE;nPAdx 


(4.49) 


the  energy  how  rate  of  each  species,  and 

-  f 

Qion  —  / 

J  A 

the  ionization  losses.  Grouping  in  Eq.  (4.47)  the  non-dominant  contributions  into  Qr,  and 
calling  W'ioo  =  the  energy  balance  becomes 

IdYd  =  W'ioo  +  Qion  +  Qr-  (4.50) 

From  the  efficiency  dehnitions  and  for  rju,r]e  ~  1,  one  has 

Qion  T  Qr 


id  -  1-1-77 


IdYd 


(4.51) 


This  expression  together  with  Eq.  (4.49)  complements  the  scaling  law  (4.34).  For  the  case  of 
Fig.  4.2,  one  has  Qt0n/ 1 dW  —  9.3%  and  Qr/ IdYd  —  3.7%  (with  qxeoo  the  main  contribution  to 
Qr)- 
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4.6  THRUSTER  PERFORMANCE 


4.6.1  Domain  of  stationary  solutions 

Model  I  showed  that  the  parametric  domain  of  stationary  solutions  was  restricted  to  r\iB  <  0 
(i.e.  <  0).  When  the  thruster  parameters  cross  the  line  riiB  =  0  the  plasma  is  unable 

to  provide  the  ion  back-flow  needed  to  keep  plasma  quasineutrality  near  the  anode,  and  an 
oscillatory  plasma  response  is  expected.  Indeed,  both  in  Ref.  [10]  and  here  we  are  being  more 
restrictive  with  the  minimum  value  of  \rjis\-  The  anode  sheath  vanishes  when  4>ab  =  0,  which, 
from  Eq.  (4.25),  corresponds  to 


~ViB 
id  ~  ViB 


(4.52) 


Although  it  is  possible  to  extend  Model  II  to  treat  cases  with  no  anode  sheath,  we  will  take 
the  case  of  sheath  collapse  as  the  practical  limit  of  existence  of  stationary  solutions.  Two  are 
the  reasons:  first,  the  limit  ViB  =  0  is  very  close;  and  second,  Eq.  (4.52)  [and  Fig.  4.4]  show 
that  the  electron  velocity  when  <f )^B  =  0  is  formally  beyond  the  limit  of  validity  of  the  diffusive 
model. 

An  important  novelty  of  Model  II  is  that  heat  conduction  introduces  an  upper  restriction  to 
\Vib\-  This  comes  from  the  Bohm  condition  at  point  B.  To  have  dvxi/dx \b  >  0  requires  Gb  >  0, 
which  using  Eq.(4.21)  and  uiB  —  0  means 


QxeB  < 


5 

2 


{rie^eVxe)  B- 


(4.53) 


Were  Gb  <  0,  the  reverse  ion  velocity  would  reach  a  minimum  subsonic  value  at  an  intermediate 
point  and  the  sonic  condition  at  point  B  would  not  be  met.  Physically,  this  means  again  the 
impossibility  to  establish  a  stationary  sheath  at  the  anode,  and  a  possible  oscillatory  response 
of  the  plasma. 

Using  Eqs.  (4.25)  and  (4.26),  Model  II  restricts  the  domain  of  of  stationary  solutions  to 


0  <  e(j)AB/TeB  <  3, 

or,  in  terms  of  plasma  currents  (and  for  xenon), 

5  •  10-3  <  \IiB\/Id  <  0.114. 


(4.54) 


4.6.2  Optimum  performance 

Figure  4.6(a)  shows  the  domain  of  stationary  solutions  in  the  plane  of  input  parameters 
(Vd,Bmax).  Figures  4.6(b)-4.6(d)  depict  the  corresponding  bands  of  output  parameters.  As 
expected  from  Eq.  (4.54)  the  domain  of  solutions  is  rather  narrow.  Figure  4.6(b)  confirms  two 
known  features  of  the  specific  impulse  Isp  =  F/rn\  first,  it  is  proportional  to  Vd '  and,  second, 
it  is  practically  independent  of  the  magnetic  field  strength  (within  the  high  efficiency  regime). 
More  interesting  to  discuss  are  the  other  plots. 

Figure  4.6(a)  shows  that  Bmax  must  increase  with  Vd  to  maintain  a  stationary  solution.  This 
same  trend  is  found  experimentally  for  optimum  thruster  performance,  which  supports  that 
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optimum  performance  is  related  to  a  non-oscillatory  response.  Accepting  that  idea,  Fig.  4.6(a) 
suggests 

Bmax  oc  vd  (4.55) 

as  the  correct  law  for  optimum  performance,  instead  of  Bmax  oc  V^2  proposed  by  Kim.  The 
discrepancy  comes  from  the  behavior  of  id  —  1,  which  Kim  assumed  independent  of  Vd,  whereas, 
in  our  case,  Figs.  4.6(c)  and  4.6(d)  suggest  the  laws 

id-l~l-770c  V~1/2.  (4.56) 

This  last  behavior  can  be  justihed  from  Eq.  (4.51).  First,  for  zero  wall  losses,  one  expects 
Te  oc  Vd:  indeed,  our  computations  yield 


Te,max  -  Vd/ 5,  (4.57) 

with  Te  max  the  maximum  electron  temperature  (but  we  could  not  find  a  precise  justification  of 
the  factor  1/5).  Then,  from  Eq.  (4.49),  one  has 

Qion  +Qr  ~  Qion  OC  Vi  OC  T^2  OC  V^2 ,  (4.58) 

assuming  ne  to  be  weakly  affected  by  Vd.  Since  —  1  is  small,  we  can  take  Id  =  rmd  ~  const 
in  Eq.  (4.51),  leading  finally  to  the  scaling  law  (4.56). 

The  scaling  law  (4.55)  and  the  increase  of  r]  with  Vd  were  obtained  also  by  Blateau  et  aim 
their  simulation  model  (with  metallic  walls)  [18].  Experimental  data  by  Manzella  et  al.  [41] 
with  a  thruster  with  ceramic  walls,  fit  well  with  B  oc  V^2  but  r]  does  not  increase  clearly  with 
Vd.  Therefore,  our  conclusions  are  that  Eqs.  (4.34)  and  (4.51)  seem  to  be  the  main  relationships 
among  thruster  parameters,  whereas  the  particular  scaling  laws  for  the  evolution  of  Bmax  and 
(1  —  rj )  with  Vd  depend  on  the  thruster  operating  point  and  the  mechanism  dominating  energy 
losses. 


4.7  CONCLUSIONS 

The  inclusion  of  heat  conduction  in  our  previous  Model  I  has  produced  the  following  changes  in 
both  the  treatment  of  the  mathematical  model  and  the  stationary  solutions:  (i)  a  new  boundary 
condition  on  the  heat  flow  at  the  anode  is  needed,  (ii)  the  singular  points  are  modified,  (iii) 
the  integration  procedure  requires  more  involved  iteration  routines,  and  (iv)  a  new  restriction 
to  the  parametric  domain  of  stationary  solutions  has  arisen.  The  new  restriction  bounds  the 
maximum  ion  reverse  flow  compatible  with  a  Bohm  transition  to  the  anode  sheath,  and  depends 
very  much  on  the  new  boundary  condition.  The  lower  and  upper  bounds  on  the  ion  reverse 
flow  leave  the  parametric  domain  of  stationary  solutions  rather  narrow.  Since  both  bounds  on 
the  ion  flow  are  related  to  the  establishment  of  a  stationary  anode  sheath,  oscillatory  responses 
are  expected  outside  the  parametric  region  of  stationary  solutions. 

Model  II  recovers  an  axial  structure  of  the  discharge  consisting  mainly  of  the  diffusion 
plus  ionization  regions,  on  the  rear  part  of  the  chamber,  followed  by  the  internal  and  external 
acceleration  regions.  The  conductive  energy  transport  turns  out  to  be  as  important  as  the 
convective  one  and  is  able  to  smooth  the  temperature  profile.  Because  of  the  short  ionization 
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region,  the  ion  temperature  is  found  to  be  negligible,  and  the  transmission  of  electrostatic 
energy  into  the  monoenergetic  ion  beam  is  very  efficient. 

The  discussion  of  the  inertial  effects  in  the  electron  dynamics  and  the  electron  drift  energy 
has  yielded  some  important  conclusions.  First,  it  has  been  shown  that  the  dominant  inertial 
term  comes  the  azimuthal  electron  motion.  Second,  this  term  can  produce  an  anomalous-like 
collisionality  in  the  near-plume,  which  could  be  an  important  feature  to  understand  the  plasma 
discharge  there;  further  investigation  on  this  issue  seems  crucial.  Third,  the  fact  that  the 
electron  inertia  can  be  significant  in  a  certain  region  does  not  imply  necessarily  the  electron 
drift  energy  to  be  important  too. 

A  detailed  evaluation  of  the  thrust  has  determined  the  external  contribution  coming  from 
the  magnetic  force  between  thruster  and  external  electrons.  This  external  thrust  is  a  fraction  of 
the  ion  acceleration  in  the  plume,  which  is  transmitted  to  the  electrons  through  the  quasineutral 
electrostatic  field.  Another  fraction  of  the  ion  acceleration  goes  to  the  plume  radial  expansion 
and  does  not  produce  thrust.  Indeed,  the  external  thrust  would  be  zero  if  the  magnetic  held 
did  not  leak  outside  the  thruster  chamber.  Therefore,  for  more  accurate  computations  it  is 
very  important  to  have  an  accurate  knowledge  of  the  profile  of  the  external  magnetic  held,  the 
approximate  position  of  the  neutralization  surface,  and  the  external  electron  dynamics. 

A  relevant  novelty  of  Model  II  is  that  the  parametric  domain  of  stationary  solutions  is 
very  narrow.  This  is  the  consequence  of  upper  and  lower  restrictions  on  the  ion  back-how 
needed  to  sustain  a  steady  solution  near  the  anode.  The  domain  of  steady  discharges  in  the 
( Bmax ,  I'd)  plane  seems  to  agree  with  the  experimental  curves  for  optimum  performance.  The 
two  basic  relationships  among  governing  parameters  at  optimum  operation  have  been  derived. 
Both  scaling  laws  Bmax  oc  Vd  and  Bmax  oc  Vd  are  possible  depending  on  the  main  mechanism 
for  energy  losses,  which  determines  the  scaling  law  for  (1  —  rj)  too. 

The  high  thrust  efficiency  obtained  by  Model  II  (about  75-85%  based  in  the  anode  mass  how) 
and  the  still  high  peak  temperature,  indicate  that  energy  losses  due  to  the  plasma  interaction 
with  lateral  walls  (in  thrusters  with  long,  ceramic  chambers,  which  are  the  main  reference 
for  our  model)  are  an  effect  of  dominant  order,  which  produces  energy  losses  as  large  as  the 
ionization  ones;  this  would  agree  with  the  results  reported  by  Kim[38].  Therefore,  efforts  made 
to  include  in  a  consistent  way  this  wallinteraction  into  Model  I  [7]  will  be  extended  to  Model 
II. 
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Figure  4.1:  Sketches  of  (a)  the  Hall  thruster,  and  (b)  the  axial  model.  Ta  =  Acnavxa,  a  =  i,  e, ... 
are  axial  flows  of  particles  of  the  different  species.  Point  A  is  the  anode,  point  P  is  the 
neutralization  surface,  and  point  S  is  the  forward  sonic  point  of  the  ion  flow.  Plasma  regions: 
AB  is  the  anode  sheath,  BD  the  diffusion  region,  DH  the  ionization  layer,  HE  the  internal 
acceleration  zone,  EP  the  near-plume,  and  Poo  the  far-plume. 
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x  (mm) 

Figure  4.2:  Axial  structure  of  the  plasma  for  xenon,  L  =  25mm,  LEp  =  17.4mm,  Ac  =  45cm2, 
'friA  =  5.3mg/s,  Vd  =  298V,  TeP  =  2.5eV,  ldaB  —  -068,  and  B(x)  as  shown  in  (a).  Output 
parameters:  Id  =  4.51  A,  F  =  0.108  N,  rj  =  80.8%,  r]iB  =  —7.1%,  4>ab  —4.4  V,  xs  =  16.4  mm, 
5e  —  25.7  deg.  The  vertical,  dashed  lines  are  the  thruster  exhaust  and  the  neutralization 
surface;  the  asterisk  situates  point  S.  In  (f)  the  solid  line  corresponds  to  the  conductive  electron 
power,  Acqxe,  and  the  dashed  line  to  the  convective  one,  Ac^Tenevxe. 
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Figure  4.3:  For  the  solution  of  Fig.  4.2:  (a)  contributions  of  different  processes  to  the  total 
electron  collision  frequency  ue;  and  (b)  ionization  frequency,  and  axial  electron  diffusion 
frequency,  vd  =  vdmelmi. 


Figure  4.4:  Profiles  of  the  Hall  parameter,  uje/ve.  the  drift-to-internal  energy  ratio,  me  <P-Tr, 
and  the  convective-to-magnetic  time  ratio,  coe  1\dvge/dx\,  for  the  case  of  Fig.  4.2. 
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Figure  4.5:  Profile  of  the  ion  temperature  for  the  case  of  Fig.  4.2. 


F*(V) 


Figure  4.6:  Band  of  stationary  solutions  for  input  parameters  Vd  and  Bmax  as  in  (a)  and  the 
rest  of  input  parameters  as  in  Fig.  4.2;  Isp  =  F/ttia ■  Lines  ’a’  and  ’b’  correspond  to  the  lower 
{4>ab  =  0)  and  upper  ( Gb  =  0)  restrictions  to 
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Chapter  5 

AXIAL  MODEL  WITH  LATERAL 
WALL  EFFECTS 

5.1  INTRODUCTION 


The  walls  of  the  discharge  channel  of  the  SPT  class  of  Hall  thrusters  is  made  of  a  ceramic 
insulator.  The  interaction  of  the  plasma  with  these  dielectric  walls  affects  greatly  the  plasma 
discharge  and  the  thruster  characteristics [46].  From  simple  reasoning  several  phenomena  are 
expected.  First,  the  zero  current  condition  at  the  dielectric  wall  implies,  in  general,  the  forma¬ 
tion  of  an  electron-repelling  sheath  there.  Second,  ions  attracted  to  the  walls  are  recombined, 
thus  reducing  the  actual  propellant  utilization  and  the  thruster  efficiency.  Third,  energetic 
electrons  arriving  at  the  wall  deposit  their  energy  there,  affecting  the  average  temperature  of 
the  plasma  in  the  channel  and,  consequently,  the  ionization  process.  Four,  ceramic  materials 
present  large  yields  for  secondary  electron  emission.  On  the  one  hand,  these  electrons  modify 
the  sheath  structure,  which  could  attain  charge  saturation  conditions [36].  On  the  other  hand, 
since  these  electrons  are  replacing  primary  electrons,  which  were  magnetically-guided,  there  is 
an  extra  near-wall  collisionality  for  the  average  electron  population[21]. 

The  effects  of  the  radial  interaction  can  be  taken  into  account  in  a  one-dimensional(lD), 
radially-averaged  model  of  the  discharge  through  appropriate  source  terms  for  the  different 
radial  phenomena  [26,  11].  However,  an  accurate  form  of  these  terms  requires  to  know  in  detail 
the  radial  plasma  response  and  this  is  not  a  simple  problem.  A  variety  of  phenomena  come 
into  play,  like  the  formation  of  a  radial  presheath/sheath  structure,  the  wall  secondary  electron 
emission  (SEE),  the  trapping  of  the  SEE  into  the  main  plasma,  the  sheath  charge-saturation, 
the  annular  geometry  of  the  chamber,  and  the  effects  of  the  axial  flow  on  the  radial  dynamics. 
References  [3,  2,  12]  present  the  development  of  a  detailed  radial  presheath/sheath  model,  which 
tries  to  include  all  this  phenomena. 

This  chapter  analyzes  solutions  resulting  from  the  inclusion  of  radial-interaction  terms  into 
the  axial  model  with  heat  conduction  of  Ch.  4.  The  modeling  of  radial  effects  is  the  same 
presented  in  Ref.  [7],  which  is  valid  for  total-trapping  of  SEE.  Related  work  from  other  authors 
is  found  in  Refs.  [25,  13]. 
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5.2  MODEL  FORMULATION 


Figure  5.1  presents  a  geometrical  sketch  of  the  axial  model;  for  the  radial  model  see  Ref.  [12]. 
We  summarize  here  the  general  hypotheses  of  the  ID  axial  model,  already  discussed  in  previous 
works  [7,  5].  In  the  ID  axial  approximation,  plasma  variables  represent  average  values  on  each 
cross-section  and  depend  only  on  the  axial  variable  x.  The  macroscopic  formulation  consists 
of  fluid-like  equations  for  ions(i),  electrons(e)  and  neutrals(n).  The  annular  channel  is  of 
length  Lc  =  xe  —  %a,  radial  area  Ac,  and  radial  width  hc.  The  plasma  discharge  extends 
outside  the  thruster  into  a  plasma  plume  of  varying  radial  area  A(x).  The  magnetic  field  is 
considered  purely  radial,  B  =  B(x)lr.  with  the  maximum  field,  Bmax,  placed  near  the  thruster 
exit,  and  different  semi-Gaussian  axial  profiles  are  used  for  B{x )  inside  the  thruster  and  in 
the  near-plume.  Electrons  are  injected  into  the  plume  at  a  neutralization  surface  (point  P) 
placed  at  a  distance  Lcat  =  xP  —  xE  from  the  channel  exit  (point  E).  The  voltage  difference 
between  anode  (point  A)  and  point  P  is  the  discharge  voltage,  14,  and  the  electron  current 
delivered  at  the  neutralization  surface  is  the  discharge  current,  Id  =  eTd.  One  part  of  this 
current  flows  outwards  and  neutralizes  the  ion  beam  current  there,  Iip\  the  other  part,  Id  —  lip, 
diffuses  inwards  across  the  radial  lines  of  the  applied  magnetic  field  and  ionizes  the  mass  flow 
of  neutrals,  friA  =  m,rra,  injected  at  the  anode.  Only  single  ionization  is  considered.  Subscript 
oo  will  refer  to  downstream  conditions  far  away  from  the  cathode  P. 

The  plasma  is  considered  quasineutral  everywhere  except  in  a  thin,  electron-repelling  sheath 
attached  to  the  anode  (region  AB  in  Fig.  5.1,  with  xE  —  xa  =  0  in  the  quasineutral  scale). 
The  macroscopic  equations  for  the  quasineutral  plasma  are  [7] 
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(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 


For  a  complete  listing  of  symbols  and  definitions  see  Ref.  [10]  and  Ch.  4.  The  ionization 
frequency  follows  =  nni?,(Te);  the  axial  diffusion  frequency  vd  for  the  magnetized  electrons 
satisfies  Eq.  (4.18)  with 

V e  T  ^ei  T  ^ano^e  T  A/; rn 

the  electron  collision  frequency,  grouping  contributions  from  e  —  n  and  e  —  i  collisions,  Bohm 
anomalous  diffusion  [35],  and  effective  wall  collisions;  vw  and  vwe  are  frequencies  accounting 
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for  particle  and  energy  losses,  to  be  defined  below;  vnw  is  the  axial  velocity  of  ions  after  wall 
recombination  and  accommodation.  TQ  =  Anava,a  =  i,e,n  are  particle  flows,  and  follow 
Tn  +  Tj  =  const  =  rm,  and  Tj  +  Te  =  const  =  IV  The  area  A{x )  and  thickness  h(x)  of  the 
plasma  jet  cross-section,  are  constant  within  the  channel  and  satisfy  Eq.  (4.19)  in  the  plume. 

The  radial  plasma- wall  interaction  is  modeled  through  the  effective  frequencies  [7]: 

—  Vwhc  \/'I \  / 771  i ,  (5.7) 

£ 

V-wm  =  fimipw)  =  Z  ?  j  (5-8) 

1  Uw 

1.65 

Vwe  -  Pe^w,  Pe(Sw)  ~  5.62  +  '  (5.9) 

1  ow 

Here,  Dw  (x)  is  a  dimensionless  parameter,  which  depends  on  the  characteristics  of  the  radial 
response  in  each  axial  position;  since  there  is  not  yet  an  accurate  way  to  estimate  it,  we  just 
take  vw(x )  constant  and  equal  to  a  value  fitting  the  experimental  data.  Functions  [im (5W )  and 
j3e(Sw )  depend  on  the  effective  secondary  emission  yield, 


MTe) 


|  Vtjtu 


Te  <  T* 
Te  >  T* 


(5.10) 


with  Tj  the  electron  temperature  of  a  quasi-Maxwellian  electron  population  for  100%  SEE, 
and  5*w  ~  0.983,  T '*  ~  0.967  T1;  the  values  of  the  SEE  yield  and  the  electron  temperature  for 
the  lateral  sheaths  to  be  charge-saturated  (yielding  [i*n  ~  60  and  j3*  ~  105).  The  value  of  Tj 
is  estimated  from  T\  ~  0 .66£'i,  with  E\  the  energy  of  a  monoenergetic  beam  yielding  100% 
SEE;  Ei  depends  on  the  wall  material  and  is  the  quantity  generally  found  in  empirical  data. 
Notice  that,  since  Te  changes  along  the  channel,  the  sheath  properties  (potential  drop,  SEE 
yield,  etcetera)  are  local  too.  In  particular,  the  lateral  sheaths  are  not  expected  to  be  charge- 
saturated  along  the  whole  channel.  For  the  axial  velocity  of  neutrals  created  by  recombination 
at  lateral  walls,  vnw .  we  assume  it  depends  on  a  wall  accommodation  factor  aw.  through  vnw  = 
T  (1  Q>w^Vxi- 


5.2.1  Boundary  conditions  and  structure  of  the  discharge 

Since  the  radial-source  terms  do  not  modify  the  mathematical  characteristics  of  the  set  of 
equations,  the  treatment  of  the  present  model  is  identical  to  the  no  wall-losses  case  developed 
in  Ch.  4.  The  eight  boundary  conditions  (5  at  point  B,  2  at  point  P,  and  one  at  the  sonic  point 
S)  and  the  numerical  integration  algorithm  are  the  same  now.  We  remind  that  these  boundary 
conditions  correspond  to  solutions  including  an  anode  sheath  AB  and  a  region  BD  of  reverse 
ion  flow  (Fig.  5.1).  This  type  of  solutions  would  represent  conditions  of  optimum  thruster 
operation.  According  to  Ch.  4,  these  solutions  exists  only  in  a  narrow  interval  of  reverse  flows, 
from  —  rjiB  ~  1%  to  —  rjiB  ~  12  —  18%,  roughly,  where  rjt  =  Yi(x)/Vm  is  the  ion  flux  ratio. 

5.3  STRUCTURE  OF  THE  DISCHARGE 

Figures  5.2  to  5.5  show  axial  profiles  of  the  main  plasma  variables  for  a  particular  operation 
point  of  a  SPT-100  type  of  thruster.  Design  and  control  parameters  are  given  in  the  caption  of 
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Fig.  5.2  and  the  profile  of  the  radial  magnetic  held  is  in  Fig.  5.2(a);  the  thruster  exhaust  is  at 
xE  =  25  mm  and  the  neutralization  surface  at  xp  =  33.5  mm. 

Points  D  and  S,  delimiting  the  ionization  region  are  placed  at  xE  ^  2.6mm  and  x$  — 
12.9  mm.  The  ion  flux  ratios  at  different  locations  are  rjtB  ~  —3%,  rys  —  79.9%,  r]iE  —  82.3%, 
rju  =  rjioo  —  83.3%.  Figure  5.3(a)  shows  the  axial  evolution  of  both  the  particle  production  and 
recombination  terms, 

Si  =  AeneVi,  Sw  =  Aenevw ,  (5-11) 

respectively;  see  and  uw  in  Fig.  5.4(a)  too.  The  total  electron  and  ion  current  lost  into  the 
walls,  Iw  =  fA  eSwdx,  amounts  to  0.77A,  whereas  the  discharge  current  is  Id  —  4.83A,  and 
the  ion  beam  current  is  7j0 0  ~  2.9A.  An  interesting  result  is  that  the  discharge  self-adjusts  the 
ionization  fraction  at  the  exit  of  the  ionization  layer  to  a  value  such  that  ion  production  and 
recombination  tend  to  be  balanced  in  the  acceleration  region,  ~  i>w.  This  implies  a  linear 
relation  between  the  wall-loss  factor  Dw  (=  0.191,  here)  and  the  ionization  inefficiency,  1  —  rju: 


V w 

1  -  Vu 


8 nrij  vnA 

7T7Tle  VnE 


hc^iE^nA") 


(5.12) 


with  Oi{Te )  =  Ri (Te) yj irme/8Te,  For  the  example  of  Fig.  5.2  the  right-hand  side  is  about  one; 
for  the  case  with  uw  =  0  discussed  in  Ch.  4  we  had  r)u  —  100%. 

Characteristic  values  of  the  electron  temperature  in  Fig.  5.2(d)  are  TeB  ~  3.4  eV,  TeD  ~ 
5.3  eV,  TeS  ~  21.9  eV,  TeE  ~  38.7  eV,  and  TeP  ~  4.4  eV.  Since  we  used  Tj  ~  39.9eV  to 
characterize  the  SEE  yield  of  the  wall  material,  energy  losses  to  the  walls  are  clearly  bounding 
the  maximum  electron  temperature  in  the  chamber.  The  axial  evolution  of  the  SEE  yield,  5W 
in  Fig.  5.3(b),  shows  that  the  two  lateral  sheaths  are  charge  saturated  (5W  —  0.983)  only  in 
a  very  tiny  region  (less  than  1mm)  around  the  chamber  exhaust.  Figure  5.3(c)  compares  the 
wall  energy  losses,  Qw  =  uweneTeAc ,  and  the  collisional  heating,  Qd  =  Udnemev^.eA.  The  large 
energy  losses  for  charge-saturated  conditions  reduce  quickly  the  electron  temperature  below  the 
CSL  value  to  a  level  where  they  can  be  compensated  by  the  Joule  heating.  This  rough  balance 
is  kept  for  most  of  the  acceleration  region. 

Figure  5.4(a)  compares  the  four  frequencies  involved  in  the  plasma  equations.  Notice  that, 
in  the  region  near  the  exhaust,  the  relation  between  uwe  and  =  Udme/rrii  reproduces  the 
relation  between  Qw  and  Qd-  Figure  5.4(b)  shows  that  the  wall  collisionality,  with  a  maximum 
of  5  x  106  s_1  at  the  chamber  exit,  is  clearly  insufficient  to  justify  the  effective  axial  collisionality 
of  the  electrons,  ne.  This  leaves  Bohm-like  diffusion  as  the  most  likely  phenomenon  explaining 
the  measured  electron  diffusion  [35]. 

Taking  the  origin  of  potentials  at  the  anode  in  Fig.  5.2(b),  potential  drops  in  the  different 


regions  are  4>ab  —  4.75  V,  4>bd  —  1-92  V,  (f>DS  —  —14.3  V,  (f>sE  —  —168.0  V,  (pEP  ~  —126.8  V, 
and  the  total  discharge  voltage  is  VJ  =  302.5  V.  Thus,  the  potential  drop  is  of  the  order  of  TeB 
in  the  sheath,  negligible  in  the  reverse  flow  region,  and  of  the  order  of  TeS/ 2  in  the  ionization 
region.  Indeed,  the  equation  for  the  ion  mechanical  energy  integrated  across  the  ionization 
region  yields 
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ViTTli{vxi  ~  Vn)dx 


(5.13) 


so  that  the  part  of  the  electrostatic  energy  not  transformed  in  beam  energy  is  3.38eV  per  ion. 
Most  of  the  potential  drop  takes  place  in  the  acceleration  region;  in  particular,  the  drop  is 
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55.6%  of  Vd  in  the  internal  region  SE  and  41.9%  in  the  external  region  EP.  In  the  near  plume, 
one  should  find  TeE  of  the  order  of  2e<f)Ep/§1  the  difference  between  these  two  magnitudes  being 
due  mainly  to  heat  conduction  effects;  in  the  present  case  we  have  2e(f>Ep/5  ~50.74eV,  so  that 
1/4  of  the  external  Joule  heating  is  transmitted  through  conduction. 

In  Ch.  4  we  reviewed  the  closed-drift,  collisional  model  for  the  electron  dynamics,  and  we 
evaluated  the  importance  of  the  inertial  effects  and  the  electron  drift  energy.  A  summary  of 
the  pertinent  electron  equations  is  presented  in  Appendix  A.  The  closed-drift,  collisional  model 
applies  if  conditions  (4.3)  and  (4.6)  apply.  Figure  5.5  plots  the  inverse  Hall  parameter,  the 
magnetic-to-convective  time  ratio,  and  the  drift-to-thermal  energy  ratio.  We  observe  that  the 
closed-drift  approximation,  Eq.  (4.3),  is  verified  in  the  whole  domain,  even  near  the  cathode, 
where  the  magnetic  held  is  vanishing.  However,  the  azimuthal  convection  and  energy  tends 
to  be  large  in  the  vicinities  of  anode  and  cathode;  a  similar  behavior  was  found  for  vw  =  0. 
Further  studies  are  needed  to  elucidate  the  global  relevance  of  these  effects. 

The  ion  temperature  should  be  of  the  order  of  electrostatic  energy  not  transformed  into  beam 
kinetic  energy.  A  computation  of  Tt  using  the  same  equation  than  in  Ch.  4  yields  TiS  ~  1  eV 
and  Tioo  ~  3.9  eV.  This  confirms  that  ion  pressure  effects  can  be  dropped  from  the  main  set 
of  equations  modeling  the  discharge.  The  comparison  of  the  present  values  of  Tt  with  those 
obtained  for  no  wall  losses  (we  had  Tioo  ~  0.11  eV  in  Ch.  4  for  a  similar  case)  demonstrates 
that  wall  losses,  by  extending  ionization  into  the  acceleration  region,  are  the  main  phenomena 
contributing  to  7). 

5.4  THRUST  AND  EFFICIENCY 

The  electrostatic  interaction  between  ions  and  electrons  plus  the  magnetic  forces  on  the  electrons 
are  the  mechanisms  by  which  plasma  thrust  is  transferred  to  the  solids  parts  of  the  engine.  Since 
the  magnetic  held  extends  some  distance  outside  the  chamber  and  the  plasma  plume  presents 
some  radial  divergence,  neither  the  momentum  thrust  of  the  plasma  at  the  chamber  exhaust 
(point  E)  nor  in  the  far  plume  (point  oc)  evaluate  correctly  the  thrust. 

Calling 

Fp(x)  =  (rriinivli  +  minnvln  +  neT^j  A(x)  (5.14) 

the  total  axial  momentum  of  the  plasma  at  a  given  axial  section, 

rP‘2 

FmPlp2  =  /  enevdeBA(x)dx,  (5.15) 

JPi 

the  magnetic  force  on  the  electron  how  between  two  channel  sections  (Pi  and  P2),  and 

poo 

F>piu  =  /  neTedA  (5.16) 

J  E 

the  axial  pressure  component  due  to  the  plume  lateral  expansion,  the  thrust  F  satishes  the  two 
relations 

F  =  FpE  +  FmEoo  =  Fp  oc  —  Dpiu.  (5-17) 
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The  leakage  of  the  magnetic  field  outside  the  thruster  channel  yields  a  contribution  to  the 
thrust,  FmEoo ,  which  is  a  unique  characteristic  of  this  type  of  electric  thrusters.  Also,  the 
balance 

Fpoo  FpE  FmEoo  T  Dpiu  (5.18) 

shows  that  one  part  of  the  external  acceleration  of  the  ions  comes  from  the  plume  lateral 
expansion  and  does  not  contribute  to  thrust.  Thus,  were  the  magnetic  field  zero  outside  the 
thruster,  all  external  acceleration  would  be  the  reaction  to  plume  radial  expansion.  For  the  case 
of  Fig.  5.2  we  have  F  ~  84.9  mN,  FpE  —  72.5  mN,  and  Fpoo  ~  86.1  mN,  so  that  the  external 
magnetic  field  contributes  a  14.6%  to  the  thrust  whereas  Dpiu  is  only  a  9%  of  the  ion  external 
acceleration. 

The  energy  balance  for  the  whole  plasma  can  be  expressed  as 


Pd  P%oo  T  Pyj  T  Pion  T  Prest  j 


where 


TioomiV2xioo/2, 


is  the  external  ion  energy  flow, 


L 


E 

Ac  I  vw 


-mi(v2xi  -  v2nw)  +  /3eTe 


nPdx 


(5.19) 


is  the  energy  lost  at  the  lateral  walls, 
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is  the  ionization  loss,  and  Prest  groups  other  contributions  (electron  heat  flow  at  A  and  oo, 
anode  losses,  energy  flow  of  neutrals...).  For  the  case  of  Fig.  5.2,  we  have  Pd  =  IdVd  —  1461  W, 
Pioo  ~  849  W,  Pw  ~  377  W,  Pion  ~  135  W,  and  Prest  —  100  W.  Thus,  P^  amounts  to  a  58%  of 
Pd ,  and  wall  losses  turn  to  be  2.5  times  larger  than  ionization  losses.  The  electron  energy  flow 
represents  90%  of  Prest- 

The  thrust  efficiency  satisfies  Eqs.  (4.44)- (4.45).  The  current  utilization  is  related  to  the 
energy  losses  through 

Pioo  Pion  T  Pw  T  Prest  /r  nn\ 

VdVe  =  -5-  =  1 - p - •  (5-20) 

1d  1d 

Since  r)e  is  very  close  to  one,  the  inefficiency  in  current  utilization  comes  from  the  energy  losses. 
For  the  case  of  Fig.  5.2,  one  has  r]  —  51.6%,  rju  —  83.3%,  r]e  —  96.9%,  rjd  —  60%,  rj0  —  106.5%. 
For  aw  =  0,  the  neutrals  have  a  non-negligible  contribution  to  the  thrust,  which  explains  that 
r]0  >  100%.  These  values  must  be  compared  with  rj  ~  81.1%,  rju  ~  100%,  and  rjd  —  85.6%  of 
the  similar  case  with  vw  =  0,  discussed  in  Ch.  4. 
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5.5  PARAMETER  INVESTIGATION 


Following  Ch.  4,  scaling  laws  involving  the  length  of  the  ionization  and  acceleration  regions, 
Las  and  LSp,  respectively,  are 


y/ aiEi/rrii 
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(5.21) 

(5.22) 


with  ud  =  Vdme/nrii  and  LAs  +  LSP  =  Lc  +  Lcat. 

Figures  5.6  and  5.7  show  the  influence  of  the  discharge  voltage  and  the  magnetic  field 
strength,  on  the  thruster  performance.  Some  plotted  parameters  are 


—  Pw/Pd,  £ w  —  Iw/Id !,  ^ ion  —  Pion/ Pd- 


(5.23) 


The  increase  of  Vd  implies  mainly  the  increase  of  the  electron  temperature,  Fig.  5.6(d).  This  is 
more  evident  outside  the  channel  where  there  are  no  wall  losses  limiting  Te.  For  Vd  >  300V,  there 
is  charge-saturation  of  the  lateral  sheaths  in  a  region  near  the  exit,  whose  extension  increases 
with  TeE  —  T\.  This  explains  that  Pw  =  twPd  increases  faster  than  Pd  and  the  decrease  of  rid- 
The  moderate  increase  of  TeS  explains  the  decrease  of  Las ,  Eq.  (5.21),  and  the  no  variation  of 
eion  =  Pion/Pd-  Since  r]u  depends  weakly  on  Te  [according  to  Eq.  (5.12)  and  within  the  range 
of  Fig.  5.6]  r)u  does  not  change  with  Vd- 

Since  Bohm-like  diffusion  is  dominant,  the  magnetic  field  strength  affects  the  axial  diffusion 
through  Ud  ~  we<Vmo  oc  Bmax.  As  Bmax  increases,  Las  in  Eq.  (5.21)  decreases;  then,  Lgp 
increases  and  r/d  in  Eq.  (5.22)  increases.  A  larger  Lgp  means  a  smaller  electric  field,  and  then 
a  lower  TeE. 

Figure  5.8  shows  the  ’optimum’  operation  line  for  the  parameters  (Vd,  Bmax).  As  expected, 
Bmax  increases  with  Vd ,  but  at  a  path  slower  than  in  the  zero-wall  case.  The  efficiency  increases 
slowly  with  Vd-  Further  analysis  of  these  results  is  needed. 

Finally,  Appendix  C  presents  studies  on  other  parameters  of  the  model.  These  numerical 
results  have  not  been  analyzed  carefully  too. 


5.6  CONCLUSIONS 

The  influence  of  different  source  terms  of  the  radial  interaction  on  the  axial  plasma  discharge 
behavior  has  been  studied.  The  solutions  obtained  with  the  resulting  model  show  a  reasonable 
agreement  with  trends  observed  in  experimental  data. 

The  SEE  yield  function  of  the  ceramic  insulator  is  crucial  for  the  energy  balance.  For  Vd 
high  enough,  the  external  heating  of  the  electron  population  leads  to  100%  SEE  in  a  region 
around  the  channel  exit.  This  means  a  very  large  energy  loss  and  a  reduction  of  the  electron 
temperature  to  a  level  where  energy  losses  at  the  walls  become  comparable  to  the  Joule  heating. 
Wall  losses  and  ionization  losses  are  the  phenomena  reducing  the  current  utilization. 

Wall  collisionality  depends  on  the  level  of  SEE  too,  but  even  for  100%  SEE,  this  collisionality 
turns  out  to  be  much  smaller  than  the  effective  collisionality  dictated  by  experimental  mea¬ 
surements.  This  leaves  Bohm-like  diffusion  as  the  main  phenomenon  explaining  the  enhanced 
electron  diffusion. 
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Since  the  discharge  presents  a  very  high  beam-energy  efficiency,  thrust  efficiency  is  deter¬ 
mined  mainly  by  the  product  of  the  current  and  propellant  utilizations.  The  current  efficiency 
is  related  directly  to  the  electron  energy  losses,  and  these  are  dominated  by  losses  at  lateral 
walls.  The  propellant  utilization  is  determined  by  the  balance  between  volume  production  and 
wall  recombination  of  plasma  in  the  acceleration  region. 

To  obtain  propellant  utilizations  in  the  range  80-90%,  the  parameter  measuring  wall  recom¬ 
bination  corresponds  to  a  value  below  the  one  corresponding  to  free,  radial  acceleration  of  the 
ions.  The  present  understanding  of  the  radial  plasma/wall  interaction  [3,  2,  12]  indicates  that 
the  only  phenomenon  able  to  inhibit  the  radial  ion  flux  would  be  a  strong  coupling  with  the 
axial  motion,  in  other  words,  the  two-dimensional  character  of  the  discharge.  Therefore,  re¬ 
search  on  the  radial  model  must  be  pursued  in  order  to  understand  better  the  mutual  coupling 
between  the  axial  and  radial  dynamics. 
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Figure  5.1:  Sketches  of  (a)  the  Hall  thruster,  and  (b)  the  axial  model.  Ta  =  Acnavxa,  a  =  i,  e, ... 
are  axial  flows  of  particles  of  the  different  species.  The  anode  is  at  point  A,  the  thruster 
exhaust  at  point  E,  and  the  neutralization  surface  at  point  P.  Region  AB  is  the  anode  sheath, 
BD  the  reverse-flow  region,  DS  the  ionization  layer,  SE  the  internal  acceleration  region,  EP 
the  near-plume,  and  Poo  the  far-plume.  Points  B,  D,  and  S  are  defined  by  ion  velocities 
vxw  =  -\jTeB/mt.  vxW  =  0,  and  vxiS  =  sjTeS/rnt. 
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Figure  5.3:  Evaluation  of  wall  losses  for  the  solution  of  Fig.  5.2.  In  last  subplot:  ’a’  is  Qw  and 
’b’  is  Qd. 


Figure  5.4:  Evaluation  of  different  collision  frequencies  for  the  case  of  Fig.  5.2. 
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Figure  5.5:  Evaluation  of  parameters  related  to  electron  dynamics  for  the  case  of  Fig.  5.2. 
Maximum  values  of  mevle/2Te  are  11.5  at  point  B  and  0.7  near  the  cathode. 
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Figure  5.7:  Influence  of  the  magnetic  field  strength  on  thruster  performance.  Rest  of  parameters 
as  in  Fig.  5.2. 
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Chapter  6 

MODEL  WITH  AN  INTERMEDIATE 
ELECTRODE 

6.1  INTRODUCTION 

Use  of  internal  electrodes  to  control  the  electric  field  and  improve  the  thruster  performance  is 
being  experimented  actively  [47,  27].  Intermediate  electrodes  control  the  discharge  profile  by 
(i)  being  biased  to  a  given  voltage  and  (ii)  exchanging  the  necessary  electron  current  with  the 
surrounding  plasma.  Active  electrodes  can  collect  or  emit  electron  currents  depending  on  the 
operational  conditions,  while  passive  electrodes  can  only  extract  electrons  from  the  plasma. 

Hall  thrusters  without  intermediate  electrode  but  operating  efficiently  are  known  to  reach 
ionization  and  beam-energy  efficiencies  around  90%.  The  thrust  efficiency  is  penalized  mainly 
by  the  current  utilization,  that  is  by  the  ratio  between  the  downstream  ion-current  and  the 
discharge  current.  Therefore,  the  role  of  an  intermediate  electrode  should  be  to  increase  the 
current  utilization  without  worsening  the  two  other  partial  efficiencies.  This  suggests  to  place  an 
electron-emitting  electrode  in  the  acceleration  region  of  the  channel  to  deliver  a  fraction  of  the 
discharge  current.  This  chapter  investigates  this  idea  theoretically  extending  the  axial  model 
presented  in  Ch.  5  to  include  the  interaction  with  the  additional  electrode.  A  detailed  energy 
equation  for  electrons  is  employed  and  the  parameter  investigation  departs  from  stationary 
solutions  with  high  thrust  efficiency  already  without  the  intermediate  electrode. 

Fruchtman  et  al.  [29]  showed  recently  that  an  electrode  placed  exactly  at  the  sonic  transition 
of  the  ion  flow  can  enhance  the  thruster  efficiency.  Some  comments  are  worth  to  be  made.  First, 
the  study  is  restricted  to  electron-absorbing  electrodes.  Second,  in  the  examples  they  present 
the  thruster  without  the  intermediate  electrode  has  a  low  thrust  efficiency,  mainly  because  of  the 
ionization  region  occupies  the  whole  channel;  by  modifying  the  electric  field,  the  intermediate 
electrode  forces  the  ionization  region  to  occupy  only  one  part  of  the  channel,  improving  thus 
the  propellant  and  the  beam-energy  efficiencies.  Third,  the  position  of  the  ion  sonic  transition 
changes  with  the  operational  conditions  and  only  for  very  particular  values  of  the  operation 
parameters,  it  will  coincide  with  the  position  of  the  intermediate  electrode.  And  four,  they  use 
a  model  with  Te  =  const  along  the  whole  channel  (and  no  reverse- flow  region,  no  plume).  The 
temperature  profile  and  the  energetic  balance  are  known  to  be  fundamental  in  determining  the 
discharge  characteristics. 
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6.2  MODEL  FORMULATION 


Figure  6.1  presents  a  geometrical  sketch  of  the  thruster  and  the  axial  discharge.  Except  for  the 
intermediate  electrode,  the  rest  of  the  model  is  identical  to  the  one  presented  in  Ch.  5.  The 
intermediate  electrode  [point  R  in  Fig.  6.1],  placed  at  x  —  is  treated  very  similarly  to  the 
external  neutralizer  cathode,  placed  at  point  P  in  the  plume.  The  bias  voltages  between  anode 
A,  electrode  R,  and  neutralizer  P  are  known;  we  call  Vd  —  4>a  ~  4>p  and  V,int  —  4>a  —  4>r-  On 
the  contrary,  the  electrical  current  exchanged  between  the  internal  plasma  and  electrodes  at  R 
and  P,  Iint  =  I<m  and  lap,  respectively,  depend  on  the  discharge  conditions  and  are  determined 
as  part  of  the  solution.  The  total  discharge  current  reaching  the  anode  is 

Id  =  IdA  =  IdP  +  IdR] 


lint  is  positive  or  negative  depending  on  the  electrode  R  emitting  or  absorbing  electrons. 

The  electron  currents  exchanged  at  points  R  and  P  introduce  discontinuities  in  some  plasma 
magnitudes.  Except  for  these  two  discontinuity  surfaces  (and  the  anode  sheath)  the  plasma  is 
continuous  and  quasineutral,  and  satisfies  the  set  of  equations 
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Here,  ud  —  u“l/ve  and  the  four  source  terms  are 

Si  =  AneUi,  SP  =  S(x  -  xP)IdP/e, 

Sw  =  Aneuw,  SR  =  S(x  -  xR)IdR/e, 


(6.1) 

(6.2) 

(6.3) 

(6.4) 

(6.5) 

(6.6) 

(6.7) 


(6.8) 


5(x)  the  Dirac  function,  and  TR  and  TP  the  temperatures  of  the  electrons  emitted  at  electrodes 
R  and  P,  respectively.  See  previous  chapters  for  the  definitions  of  the  rest  of  symbols  and 
variables.  Notice  that  only  the  particle  and  energy  equations  of  the  electrons  are  affected  by 
the  exchange  of  current.  On  the  contrary,  for  a  collisional  electron  motion,  ud  and  the  Ohm’s 
law  are  unaffected. 
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To  integrate  the  equations  across  point  R  the  jump  conditions  there  must  be  defined.  Clearly, 
Eq.  (6.4)  states  that  the  jump  in  the  electron  flow  is 

^eR  ~  r eR  =  ^ dR/e ■  (6-9) 

The  next  issue  is  whether  a  discontinuity  in  the  plasma  potential  takes  place  at  point  R.  This 
would  mean  the  presence  of  a  double  layer  (i.e.  a  double  Debye  sheath  with  total  electri¬ 
cal  charge  equal  to  zero).  Although  a  detailed  study  of  the  Poisson  equation  has  not  been 
completed,  the  analysis  of  the  plasma  conditions  suggests  that  a  double  layer  cannot  develop 
there.  Then,  we  have  q bR  =  q and  the  ion  equations  yield  v+iR  =  v~iR  and  n+R  —  n~R.  As  a 
consequence,  the  jump  in  electron  current  implies  a  jump  in  electron  velocity: 

VxeR  ~  VxeR  =  W <™eR-  (6-10) 

Finally,  using  Eqs.  (6.6)  and  (6.9),  the  jump  of  the  heat  conduction  flux  is 

(«£*  -  9 =  |(Tb  -  T.R)!f.  (6.11) 

Here,  TeR  is  the  temperature  at  point  R  of  the  (unique)  electron  population,  whereas  TR  is  the 
temperature  of  the  electrons  exchanged  between  the  plasma  and  the  electrode.  Thus,  one  has 
Tr  =  TeR  for  an  absorbing  electrode,  whereas  TR  is  given  for  an  emissive  electrode  (we  will  take 
Tr  r-j  TP).  Notice  that  at  the  emissive  neutralizer  P  we  take  TeP  =  TP  since  all  electrons  there 
come  from  the  cathode. 

The  rest  of  boundary  conditions  needed  to  integrate  the  plasma  equations  are  identical  to 
the  case  without  electrode  R.  The  integration  procedure  depends  on  the  position  of  the  ion  sonic 
point  S  (defined  by  vxis  —  \jTeS/rrii)  relative  to  the  electrode  R.  That  position  is  determined 
as  part  of  the  solution  and  marks  the  beginning  of  the  acceleration  region.  In  this  chapter  we 
solve  cases  where  point  S  is  upstream  of  the  electrode  R.  On  the  one  hand,  they  allow  us  to 
use  the  same  integration  scheme  than  in  the  previous  chapter.  On  the  other  hand  and  as  we 
commented  above,  these  cases  should  correspond  to  the  appropriate  configuration  to  maximize 
the  thruster  efficiency. 

We  define  the  relative  flows  of  the  different  species  as 

r]a{x)  =  Ta(x)/Tm.  (6.12) 


To  compare  the  different  cases  we  define  the  thrust  efficiency  as 


F2 

2m a  Pd 


VeVuVdVo  flint  i 


with  F  the  thrust,  defined  as  in  Ch.  5, 


(6.13) 


Pd  —  IdpVd  T  IdRYint  —  IdYd  lintiYd 


(6.14) 


the  power  spent  in  the  discharge, 


flint  —  IdVd/Pd 


(6.15) 
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the  power  gain  due  to  the  intermediate  electrode,  and 


Vu  = 


mA 


Vd  = 


c  ±  too 


h 


Ve  = 


rriiV 


2eVd 


Vo  = 


{TimiVxiyj 


(6.16) 


the  propellant  utilization,  the  current  utilization,  the  beam-energy  efficiency,  the  other-effects 
efficiency  (plume  radial  expansion,  thrust  of  neutrals,  anode  sheath,...),  respectively. 


6.3  STRUCTURE  OF  THE  DISCHARGE 

Figure  6.2  shows  a  solution  for  a  conventional  mid-power  thruster  without  intermediate  elec¬ 
trode  and  for  near-optimum  operational  conditions.  Then,  an  intermediate  electrode  is  placed 
at  xr  —  18  mm,  and  solutions  are  sought  departing  from  the  one  shown  in  Fig.  6.2  (which  has 
lint  —  0  and  Vint  —  —  4>r  —  52  V).  Figure  6.3  shows  the  solution  for  Vint  —  73.5  V,  when 

the  electrode  R  is  emitting  Iint  —  0.5  A.  Comparing  the  thruster  performances  for  the  cases  of 
Figs.  6.2  and  6.3,  we  find  a  relative  gain  of  3%  in  discharge  power  and  of  5%  in  efficiency,  when 
a  10%  of  the  discharge  current  is  emitted  by  electrode  R.  Thrust  and  propellant  utilization  do 
not  change  practically. 

Notice  in  Fig.  6.3  that  only  re  and  dTe/dx  are  discontinuous  at  R.  By  decreasing  <j)R  we  are 
forcing  lower  potentials  and  electric  fields  in  the  plume;  the  electron  current  between  R  and 
P  decreases  with  Ilnt  too.  The  consequence  is  a  smaller  electron  temperature  at  the  thruster 
exhaust  with  the  subsequent  reduction  of  the  wall  energy  losses.  We  showed  in  Ch.  5  that  these 
losses  are  related  directly  to  the  current  inefficiency,  1  —  Iioo/Idp ;  from  Fig.  6.2  to  Fig.  6.3, 
hoc/IdP  has  increased  a  5%.  This  increment  explains  the  improvement  of  rj ,  since  the  rest  of 
partial  efficiencies  change  very  little;  in  particular,  the  variation  of  rju  is  small  because  of  the 
electrode  R  is  downstream  of  the  ionization  region. 


6.4  PARAMETER  INVESTIGATION 

Steady  solutions  with  a  reverse-flow  region  exist  only  in  the  narrow  band  1%  <  1 77^^ |  <  15%, 
roughly(Ch.  5).  Figures  6.2  and  6.3  show  the  increase  of  the  reverse  flow  with  Iint.  That  narrow 
parameter  band  would  represent  ’optimum’  points  of  operation,  as  we  discussed  in  Ch.  5.  To 
keep  the  thruster  operating  in  that  band,  we  look  for  pairs  of  control  parameters  yielding- 
solutions  which  keep  rjiB  ~  const  within  that  range.  Figure  6.4  shows  the  optimum  operational 
points  of  the  pair  {Vint,  Bmax ),  —77,#  =  2.5%,  and  the  rest  of  parameters  as  in  Fig.  6.2;  notice  that 
lint ,  although  represented  in  the  abscissa,  is  an  output  parameter.  Both  absorbing  ( Iint  <  0) 
and  emitting  ( Iint  >  0)  electrodes  have  been  considered;  Iint  —  0  correspond  to  the  case  of  Fig. 
6.2.  The  thruster  efficiency  increases  with  the  electron  emission  at  R,  provided  that  magnetic 
field  is  increased  too;  7/  goes  from  51.5%  to  60%,  for  Iint  changing  from  0  to  1.3  A  and  Bmax 
from  233  gauss  to  362  gauss.  Since  rj/f3int  and  the  total  discharge  current  IdA  vary  very  little 
with  Iint,  the  gain  in  efficiency  comes  directly  from  the  gain  in  discharge  power, 

Let  us  try  to  understand  now  the  existence  of  a  maximum  of  r\  (for  Iint  ~  1.3  A).  The 
lower,  ideal  limit  for  the  electron  current  delivered  by  the  cathode  is  IdP  —  Up-  Since  IdA 


92 


and  lip  remain  almost  constant  with  Iint,  the  difference  I^a  ~  Up  yields  the  maximum  value 
of  the  intermediate  current.  Figure  6.4  shows  that  the  derivative  —dVint/dIint  increases  as  Iint 
approaches  that  maximum.  The  consequence  is  that  Iint{Yd~  Vint )  in  Eq.(6.14),  and,  therefore, 
flint  and  T)  reach  a  maximum. 

For  the  parameter  evolution  of  Fig.  6.4,  point  S  moves  downstream  with  Iint  increasing;  for 
lint  —  1-2  A  it  is  placed  at  xs  ~  10  mm.  According  to  the  law 

_  r  1  /cTj-Fj  V{ 

%s  =  Las  ~  — \ - 

V  me  Vi 

the  displacement  of  point  S  is  due  to  the  change  of  Bmax. 

Figure  6.4  shows  that  the  case  Ijnt  <  0,  corresponding  to  an  electron-absorbing  electrode,  is 
clearly  disadvantageous,  at  least  when  the  thruster  is  operating  efficiently  without  intermediate 
electrode.  As  Iint  becomes  more  negative,  point  S  approaches  point  R.  The  situation  considered 
by  Fruchtman  et  ai,  with  the  sonic  point  just  at  the  electrode  R  is  reached  for  a  particular 
value  of  the  intermediate  current  around  -1.2  A;  the  thrust  efficiency  then  is  43%  only. 

Figure  6.5  shows  the  structure  of  the  discharge  for  the  case  of  maximum  efficiency  in  Fig. 
6.4.  Here,  the  phenomena  we  observed  in  Fig.  6.3,  are  more  pronounced.  The  most  evident 
one  is  the  large  modification  of  the  temperature  profile.  The  low  Joule  heating  between  R 
and  S  yield  low  electron  temperatures  there.  For  this  case,  lateral  sheaths  are  far  from  being 
charge-saturated.  Indeed,  the  increment  of  the  electric  field  between  R  and  S,  yield  a  large 
Joule  heating  which  produces  the  intermediate  maximum  of  Te,  but  this  one  is  still  below  the 
charge-saturation  limit. 

6.5  CONCLUSIONS 

A  previous  model  of  the  plasma  discharge  has  been  extended  to  include  the  exchange  of  electrical 
current  with  an  internal  electrode  biased  to  a  given  voltage.  The  electrode  can  either  emit  or 
absorb  electrons.  Only  configurations  with  the  electrode  placed  in  the  plasma  acceleration 
region  have  been  treated  but  these  are  thought  to  be  the  most  interesting  ones  in  practice. 

It  has  been  demonstrated  that  the  efficiency  of  a  mid-power  thruster  can  be  increased  from 
about  50%  to  60%  with  1/4  of  the  electron  current  emitted  from  the  intermediate  electrode.  The 
increase  in  efficiency  is  the  consequence  of  a  reduction  in  discharge  power,  for  given  discharge 
voltage  and  mass  flow.  When  the  electrode  potential  is  decreased,  the  plasma  increases  the 
exchange  of  the  intermediate  current,  whereas  the  total  discharge  current,  the  thrust,  and  the 
propellant  utilization  are  kept  almost  constant.  The  consequences  are  lower  electron  heating, 
lower  wall  losses,  and  larger  utilization  of  the  neutralizer  current. 

A  parameter  investigation  on  the  pair  magnetic-field-strength/intermediate-emitted-current 
concludes  that  Bmax  must  be  increased  with  Iint  to  follow  points  of  stationary,  efficient  oper¬ 
ation.  It  also  demonstrates  that  there  exist  an  intermediate  current  which  gives  maximum 
efficiency.  This  maximum  is  related  to  the  fact  that  the  total  discharge  current  varies  weakly 
with  the  intermediate  current,  and  this  one  is  upper  bounded  by  the  cathode  current  required 
to  neutralize  the  ion  flow  in  the  plume.  Finally,  it  has  been  shown  that  an  absorbing  electrode 
placed  in  the  acceleration  region  has  only  negative  effects  on  the  thruster  performance,  when 
the  thruster  is  operated  correctly.  Further  optimization  studies  are  needed. 


(6.17) 
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Figure  6.1:  Sketches  of  the  axial  model.  The  anode  is  at  point  A,  the  thruster  exhaust  at 
point  E,  the  neutralization  surface  at  point  P,  and  the  intermediate  electrode  at  point  R.  The 
ionization  region  will  correspond  to  the  region  between  point  D  of  zero  velocity  and  the  sonic 
point  S.  Region  AB  is  the  anode  sheath.  The  voltages  Va  and  Vint  are  known,  whereas  Er  and 
IdP  are  determined  by  the  thruster/plasma  response.  Ta  =  Acnavxa,  a  —  i,  e, ...  are  axial  flows 
of  particles  of  the  different  species. 
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Figure  6.2:  Axial  structure  of  the  plasma  discharge.  Input  parameters:  Lc  =  25mm,  Lcat  = 
8.5mm,  Ac  =  40cm2,  Vd  =  296V,  m  =  4.78mg/s,  TeP  =  4.4eV,  16o:ano  ~  0.1986,  T\  ~  39.93eV, 
vw  ~  0.191,  aw  =  0.  The  intermediate  electrode  is  placed  at  xr  =  18  mm  and  is  biased  at  a 
voltage  Vint  =  §A  -  (f)R  =  52  V. 
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Figure  6.3:  Plasma  response  with  an  intermediate  electrode  at  xr  =  18  mm,  biased  to  = 
^a~  4>r  =  73.5  V,  and  emitting  Iint  =  0.5  A.  The  rest  of  parameters  are  as  in  Fig.  6.2.  (where 
<t>A  -  <t>R  =  52  V). 
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Figure  6.4:  Optimum  operation:  Values  of  Vint  and  Bmax  keeping  —rjiB  =  const  =  2.5%  and 
the  electrode  at  xr  =  18  mm.  The  rest  of  parameters  are  as  in  Fig.  6.2.  Iint  <  0  correspond  to 
an  absorbing  electrode,  I.mt  =  0  to  no  electrode,  and  Iint  >  0  to  an  emitting  electrode. 
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Appendix  A 

DERIVATION  OF  THE  RADIAL 
PRESHEATH  MODEL 


In  a  Hall  thruster,  ions  are  created  by  electron  impact  of  a  neutral  gas,  most  of  them  are 
singly-charged,  and  magnetic  effects  on  them  can  be  ignored.  From  the  whole  axisymmetric 
model,  the  equations  relevant  to  study  radial  ion  dynamics  are 
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Here  is  the  ionization  frequency,  Tn  is  the  temperature  of  the  neutral  gas,  and  the  rest 
of  symbols  are  conventional.  Ion  heat  conduction  is  difficult  to  estimate;  here  it  has  been 
disregarded,  based  in  the  conclusions  of  Scheuer  and  Emmert[50].  In  Hall  thrusters,  Tn  is 
negligible  compared  to  the  plasma  temperature  Te;  here,  it  is  kept  to  compare  our  model  with 
previous  ones. 

From  the  preceding  equations,  in  each  axial  section  of  the  chamber,  ion  radial  dynamics 
verify 
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(A.4) 

(A.5) 

(A.6) 


where  axial  dynamics  are  included  through  three  frequency  terms, 


1  d 

v,„  —  vi - —  ( riiVxi ), 


Vr  ~  Vi  + 


Hi  dx 

h xi  dvri 
Vri  dx  ’ 


Vi  (l  2  3  \  t> 

T,  l2mi^  +  2r"J-2 


Vt 


5  dvT;  3vrj  drijTj 


2  Ox  2n,T,  dx 


(A.T) 


These  effective  frequencies  are  unknown  functions  of  x  and  r.  The  ID  radial  model  is  based 
in  taking  Vj  =  Vj(r),  ( j  =  i,w,r,t ),  as  known  functions  in  a  generic  axial  position  x.  [Effective 
frequencies  for  charge-exchange,  electron-neutral,  and  electron-ion  collisions  have  been  omitted 
in  Eq.  (A. 2),  based  in  the  assumptions  vcx  <  0(vw )  and  ven ,  vei  <C  uwmi/me,  holding  generally 
inside  a  Hall  thruster.  It  would  be  no  problem  to  include  them  in  the  definition  of  ur.] 

It  is  worth  noting  that  the  equations  defining  uw  and  Eq.  (A. 7),  also  correspond  to  the 
axial  equations  of  continuity  and  energy  for  the  ’average’  (i.e.  r-independent)  plasma  state  in 
an  axial  cross-section, 


dx^d^xi)  —  {Vi  Vyj^Tli, 

5/3  x  dvxi 

^XiV^+nX— 


=  Viidnintvli  +  -  ('.n.X, 


(A.8) 

(A-9) 


Since  the  losses  of  internal  energy  at  lateral  walls  per  particle  are  5Tj/2,  vt  should  be  equal  to 
5uw/2,  which  coincides  with  the  result  found  in  Sec.  II,  Eq.  (1.12),  from  the  analysis  of  the 
radial  energy  equation. 

In  a  Hall  thruster,  electrons  diffuse  axially  in  a  quasi-closed  drift  induced  by  a  radial  mag¬ 
netic  field,  whereas  the  remain  confined  radially  by  the  self-induced  electric  potential.  There¬ 
fore,  we  can  assume  that  in  the  radial  direction  electrons  are  in  Maxwell- Boltzmann  equilibrium 
at  a  temperature  Te, 


m  one  oo 

0  ~  Te  —jr  +  ene  — 
dr  dr 


The  Poisson  equation  for  the  axisymmetric  problem  if 


(A. 10) 


a2<i, 

dx2 


1  d  /  d(j)\ 
^  r  dr  \  dr  / 


—  K  -  rii). 
eo 


(A.ll) 


Since  the  Debye  length  is  much  smaller  than  both  the  thruster  dimensions  and  the  relevant 
mean  free  paths,  the  zero  Debye-length  limit  can  be  adopted  to  solve  the  problem.  This  means 
to  consider  the  plasma  quasineutral  (ne  ~  n*)  everywhere  except  in  thin  space-charge  sheaths 
around  the  chamber  walls,  where  a  ID  form  of  Eq.  (A.ll)  can  be  used. 
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Appendix  B 

ON  THE  PRESHEATH/SHEATH 
TRANSITION 


The  quasineutral  equations  (1.1)-(1.5)  show  that  it  is  impossible  to  have  \vri\  >  cs  in  the 
presheath.  But  these  equations  by  themselves  do  not  forbid  a  presheath/sheath  transition  with 
\vri\  <  cs ,  as  Keidar  et  al. [37]  propose.  To  settle  that  | vri  =  cs  is  the  only  possible  transition 
condition  is  necessary  to  include  the  Poisson  equation  in  the  analysis. 

Let  us  consider  that  the  quasineutral  solution  extends  until  certain  point  Q  with  urjg  <  csq 
and  neQ  =  rqg.  The  possibility  of  a  transition  to  a  non-neutral  sheath  depends  on  the  local 
behavior  of  the  plasma  and  Poisson  equations  around  point  Q.  Taking  into  account  that  the 
sheath  is  a  thin  region  tied  to  the  wall,  radial  gradients  would  dominate  within  it  and  in  the 
transition  region,  even  for  a  plasma  flowing  axially.  Then,  the  Poisson  equation  (A.  11)  reduces 
there  to  the  planar  form 

7^  =  nM)  ~  n*(^r)>  (R1) 

where  the  two  density  functions  come  from  Eqs.  (1.1)-(1.5).  Calling  A (j>  =  —  A r  =  r  —  tq, 

and  XdQ  =  \J 6oTe/e2neQ1  the  local  expansion  of  Eq.  (B.l)  around  point  Q  is 


with 


+  aA<^  + b A(f>2  - cAr’ 


a  = 


miVriQ  ~  3TiQ 


-1, 


b  = 


-( 

Un  V 


d2rij 


d2ne 


2eriiQ  V  d(j)2  d(f)2 


c  =  — 


7’  dm 


Q 


eniQ  dr 


(B.2) 


(B.3) 


Equation  (B.2)  shows  that  for  a  >  0,  A^(Ar)  can  present  only  spatial  oscillations  around  the 
quasineutral  solution,  aAcj)  ~  cAr.  Therefore,  a  transition  to  a  non-neutral  sheath  is  only 
possible  for  a  <  0. 

For  cold  ions  (Tjg  =  0),  a  <  0  is  equivalent  to  vriQ  >  csq.  For  warm  ions,  one  has  a  <  0 
for  vriQ  >  csq  and  for  vriQ  <  but  in  this  last  case  the  sheath  equations  would  not 

admit  to  have  >  7 hTiQ/3mi  anywhere  and  wall  conditions  could  not  be  fulfilled  [at  least 
for  Tjg/Te  <  0(1)].  Therefore,  a  transition  to  a  sheath  structure  is  only  possible  for  |vri|  >  cs. 
To  conclude,  since  urjg  >  csq  cannot  be  reached  by  the  quasineutral  solution,  the  sonic  Bohm 
condition,  vriQ  =  csq,  (i.e.,  a  =  0)  is  the  only  possible  transition  condition. 


101 


Once  the  transition  point  is  uniquely  determined,  the  solution  of  Eq.  (B.2)  with  a  =  0 
yields  the  transition  behavior  of  the  electric  potential  and  field.  This  problem  was  solved  al¬ 
ready  by  Tonks  and  Langmuir  (see  Fig.  6  in  Ref.  [54])  and  with  some  more  detail  and  generality 
by  Lam [48];  both  works  use  kinetic  models  but  the  deduction  is  the  same  for  a  macroscopic 
approach.  The  solution  of  the  transition  region  shows  that  (i)  the  electric  field  increases  con¬ 
tinuously  and  smoothly  from  d(j)/dr  ~  Te/eh  at  the  presheath  side  to  d<f>/dr  ~  Te/eXd  at  the 
sheath  side  and  (ii) 


d(e<f>/Te) 

d( 


=  0 


(B-4) 


is  the  correct  asymptotic  boundary  condition  at  the  sheath  edge  in  the  sheath  distinguished 
scale,  (  =  (r  -  rw)j Xd. 

In  the  presence  of  a  magnetic  field  oblique  to  the  wall,  another  situation  where  Beilis  and 
Keidar  proposed  a  subsonic  ion  flow  at  the  sheath  edge[14],  Ref.  [1]  demonstrated  that  the  sonic 
Bohm  condition  continues  to  be  the  correct  presheath/sheath  transition. 
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Appendix  C 

PARAMETER  INVESTIGATION  ON 
THE  AXIAL  MODEL 


Here  we  present  graphics  showing  the  influence  on  the  plasma  and  thruster  performances  of 
several  of  the  parameters  of  the  axial  model  presented  in  Chapter  5.  In  all  figures  the  base 
case  is  that  of  Fig.  5.2  and  only  the  parameters  shown  in  the  abcissae  is  varied.  Chapter  5 
already  showed  the  influence  of  and  Bmax.  The  rest  of  cases  is  presented  here  because  they 
have  not  been  analyzed  in  detail  yet.  Also,  the  reader  must  be  aware  that  in  several  figure 
labels,  the  name  of  certain  magnitudes  is  not  the  same  we  used  in  the  main  text. 
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Figure  C.l:  Influence  of  the  anode  mass  flow;  Tm  is  actually  m  = 
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Figure  C.2:  Influence  of  the  velocity  of  neutrals  at  the  anode 
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Figure  C.7:  Influence  of  the  channel  length:  the  position  of  the  anode  is  changed,  the  B-field 
profile  remains  constant  with  respect  to  the  channel  exhaust. 


Figure  C.8:  Influence  of  the  position  of  the  neutralization  surface;  xcat  is  LEP  in  the  main  text. 
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Figure  C.9:  Influence  of  the  axial  length  scale:  all  axial  magnitudes  (Lcat,  xBmax,  LBint ,  and 
LBext)  are  scaled  proportional  to  L,  but  radial  magnitudes  are  not  modified. 


Figure  C.10:  Influence  of  the  channel  width;  h  is  hc  in  the  main  text. 
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